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A key feature of the Brown and York definition of quasilocal energy is that under local 
^ \ boosts of the fleet of observers measuring the energy, the quasilocal energy surface density 

transforms as one would expect based on the equivalence principle, namely, like E in the 
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special relativity formula: E — p = m . In this paper I will provide physical motivation 



for the general relativistic analogue of this formula, and thereby arrive at a geometrically 
natural definition of an 'invariant quasilocal energy', or IQE. In analogy with the invariant 
^ \ mass m, the IQE is invariant under local boosts of the fleet of observers on a given two- 

surface S in spacetime. A reference energy subtraction procedure is required, but in contrast 
to the Brown- York procedure, S is isometrically embedded in a /emr-dimensional reference 
^ ■ spacetime of one's choosing. For example, it is well known that any sphere, round or not, 

can always be isometrically embedded into Minkowski space, even if its scalar curvature 
is not everywhere positive. So rather than embeddability being a concern, the problem 
now is that such embeddings are not unique, leading to an ambiguity in the reference IQE. 
However, in this co dimension- two setting there are two curvatures associated with S: the 
curvature of its tangent bundle, and the curvature of its normal bundle. Taking advantage 
of this fact I will suggest a possible way to resolve the embedding ambiguity, which at the 
same time will be seen to incorporate angular momentum into the energy at the quasilocal 
level. I will analyze the IQE in the following cases: both the spatial and future null infinity 
limits of a large sphere in asymptotically flat spacetimes; a small sphere shrinking to a 
point along either spatial or null directions; and finally, in asymptotically anti-de Sitter 
spacetimes. The last case reveals a striking similarity between the reference IQE and a 
certain counterterm energy recently proposed in the context of the conjectured AdS/CFT 
correspondence. 



1 E-mail address: epp@rri.emet. in 



1 Introduction 



It is generally agreed that gravitational energy exists, but because of the equivalence 
principle it cannot be localized. The notion of quasilocal energy is currently one of the most 
promising descriptions of energy in the context of general relativity, and can be characterized 
simply as follows. The total energy, including both matter and gravitational contributions, 
contained in a finite spatial volume £ can be defined only as the integral of some energy 
surface density over its two-surface boundary, S = dH. This implies that, strictly speaking, 
there is no such thing as a local energy volume density, except that which arises from the 
small S limit of quasilocal energy^] And even this local notion is not truly local because it 
cannot be integrated over a finite volume unless one is willing to ignore effects due to gravity. 
In short, energy is associated with closed spacelike two-surfaces in spacetime, not points. 

There is also a growing consensus that the ADM and Bondi-Sachs masses are simply not 
enough. We need some definition of energy that is 'more local' than these, i.e., a quasilo- 
cal definition that does not rely on the existence of an asymptotically flat region 0. For 
example, recent proponents of this movement are Ashtekar et al [0,0, who have introduced 
the quasilocal idea of an 'isolated horizon' to describe a black hole. They articulate several 
reasons for this need, and it is useful to paraphrase here at least part of their argument: 
Let us accept that a black hole is a thermodynamic object, and so obeys the first law: 
5E — T SS + ■ ■ ■. Now suppose that the universe is asymptotically flat in spatial directions, 
and contains a single black hole. Then E in the first law is the ADM mass. But if there is 
anything else in the universe then E is not the ADM mass, and the question arises, What 
expression is to be used for E in the first law? In other words, we expect that we can put 
something else in the universe, say a galaxy somewhere, such that the black hole we started 
with, considered 'by itself, will still behave as more or less the same thermodynamic object, 
with the same mass, radiating at the same temperature as before, and with the same entropy 
equal to one quarter its area. This expectation requires the ability to compute the energy of 
a given system contained within a finite closed surface, rather than merely the total energy 
of all such systems comprising the whole universe. 

Thus quasilocal energy lies between the notions of local energy density and total energy 
of an isolated system, in the sense that it is expected to give the energy contained in any 
volume, no matter how small or large. Although the equivalence principle precludes the 
existence of a local gravitational energy density, it does not prevent us from evaluating the 
(quasilocal) gravitational energy in an arbitrarily small but nonvanishing volume S. This is 
because no matter how small £ is, S — <9£ is not a point, but rather the boundary of some 
neighborhood of a point, and so we are always inherently making a 'tidal force measurement'. 
In this sense quasilocal energy is distinct from attempts to define a local gravitational energy 
density based on certain symmetries of the action, and the concept of a Nother charge.0 At 
the other extreme is the Komar mass (or the closely related ADM mass), i.e., the total 

X A notable exception is the Tolman density, which integrates to the Komar mass [Q. But it can be 
defined only when the spacetime possesses special properties, namely a timelike Killing vector field and an 
asymptotically flat spatial infinity, and so tells us little about the nature of energy in a general context. 

2 A recent discussion of the connection between pseudotensor methods and the quasilocal idea can be found 
in Ref. §. 
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energy associated with the time translation symmetry of an isolated system. As emphasized 
in Ref. , this gravitational conserved charge is intimately connected with a lapse function, 
whereas quasilocal energy need not make any reference to a lapse function. The point is, 
the two are conceptually distinct |2],[(J , even though in some circumstances one might expect 
their numerical values to coincide. 

Currently there are several contenders for a good definition of quasilocal energy (see 
Refs. [[{], HH and the references therein). The two that interest us at the moment are 
Brown and York's 'canonical quasilocal energy' (or CQE) |7|], and the various definitions 
based on the integral over S of the Witten-Nester two-form (the two-form used in Witten's 



proof of the positive energy theorem |10[|Tl|]). The latter approach uses spinorial methods, 
and the different definitions are distinguished by the choice of supplementary equation the 
S-spinors are supposed to satisfy, for example the Sen- Witten equation |IH|l^| , the Dougan- 



Mason equation |T3| , or the Ludvigsen-Vickers equation [jTJj]. The Brown- York definition of 
quasilocal energy has the form[] 

CQE = -— / dSk-CQE rel , (1.1) 
8tt Js 

in geometrized units, with G — c — 1. The CQE is supposed to be the energy of the gravita- 
tional and matter fields contained in a finite spatial volume S, whose boundary two-surface 
is S = <9S. dS is the induced integration measure on S, and k is the trace of the extrin- 
sic curvature of S as embedded in S. Thus, —k/(8ii) is the Brown- York quasilocal energy 
surface density. When E is asymptotically flat the integral in Eq. (|1 . 1| ) (the 'unreferenced' 
CQE) diverges as S is taken to infinity, and a reference term, denoted CQE ref , is required 
to regulate the energy. Brown and York's prescription is to choose 

CQE ref = — / dSk iei , (1.2) 
8ir Js 

where k rei is the trace of the extrinsic curvature of an isometric embedding of S into some 
reference space, usually taken to be flat 1R 3 . With this choice the resulting CQE reduces 
to the ADM mass when S is taken to infinity J/J. While the CQE has a host of desirable 
properties, neatly summarized in Ref. [|15|], the embedding prescription needed to evaluate 
CQE ref is not entirely satisfactory, because not all two-surfaces that arise in practice can 
be embedded into flat 1R 3 . A ready example is the horizon of the Kerr black hole, which 
fails to be embeddable in flat 1R 3 when the angular momentum exceeds the irreducible mass 
(but is not yet an extremal black hole), and the two-sphere develops regions with negative 
scalar curvature ||16|| . While this is but one example, it is noteworthy that the breakdown 



of embeddability is in this case associated with angular momentum and negative scalar 
curvature. It is precisely such issues: embeddability, angular momentum, and negative 
scalar curvature, that will figure prominently in this paper, and will be seen to be subtly 
intertwined. 

Although a relationship between the Brown- York quasilocal energy and the spinorial 
definitions based on the Witten-Nester integral is not immediately obvious, Lau |15[ has 

3 We use a sign convention for extrinsic curvatures opposite to that of Brown and York, hence the negative 
sign in front of this integral. 
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shown that spinors may always be chosen so that the resulting spinorial definition is equal 
to the unreferenced Brown- York quasilocal energy in Eq. ( |1 . 1| ) . Moreover, he shows that the 
role of the Sen- Witten equation is to provide a definite reference point for the energy, which 
is not in general the same as CQE rcf in Eq. (JO]). The point beine; made here is two-fold: 
(i) the unreferenced Brown- York quasilocal energy seems to be robust, and (ii) all of the 
problems lie in choosing a suitable reference energy. The various prescriptions are either not 
generally well defined, or they do not agree with each other. 

I will now present a brief review of the Brown- York approach in a form that will be useful 
to us later. The classical energy-momentum tensor of matter is a local concept, associated 
with a spacetime point. It is defined for any field theory residing on a nondynamical back- 
ground spacetime (M, g) via the functional derivative of the (first order) matter action with 
respect to the metric, as follows: 



2 5 9 I™^g}= / AV^C^. (1-3) 

J M 

Here <p denotes the matter field(s) in question, and the factor of two on the left is a conven- 
tion. Usually, as we will assume here, there is no boundary term arising from this variation 
(for 'minimally coupled matter'), but in case there is it does not change the essence of the 
following argument, it just adds an interesting dimension to it. T^ at so defined is covari- 
antly conserved, as follows from the matter Euler-Lagrange equations. This prescription for 
learning about matter energy-momentum gives reasonable answers for all field theories, and 
so it is natural to try the same thing for gravity. In this case one finds, for the usual first 



order action [17], 



2 5 g I^[g] = I ~^G ab ) 5g ab + I d A x^[ -^-U ab ) 5 lah . (1.4) 
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Inspecting the bulk term one is thus tempted to define T g f av := -G ab /(% 7r) as the local 
energy-momentum tensor of the gravitational field, where G ab is the Einstein tensor. And 
this is perfectly reasonable: it is covariantly conserved — in this case identically so, via the 
contracted Bianchi identity. Moreover, its on-shell value is zero, in full accord with the 
equivalence principle, i.e., there is no nontrivial local energy- momentum tensor for the grav- 
itational field. 

In fact this absence of a nontrivial local energy-momentum tensor is true not only for 
the gravitational field, but also for any system comprised of both matter and gravity. To see 
this we need only make the spacetime metric dynamical, in which case the matter action in 
Eq. ( |1.3|) must be augmented by the gravity action. Adding Eqs. ( |1.3|) and ( |1.4j ) one finds 



for the total action 



2 5 g I tot [cp, g] = J d'x J=g (t± - 8g ab + jf d 3 x ^ (-^ ab ) 6 lab . (1.5) 



Thus one is led to identify T^ b t := T^ ai — G ab / (8ir) as the total local energy-momentum tensor 
for matter plus gravity. It has the desirable property of being covariantly conserved, but turns 
out to be just zero by the Einstein equations. If this argument is taken seriously we learn 
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that, as soon as we add gravity to any matter system, the notion of a nontrivial local energy- 
momentum tensor disappears. Furthermore, one might interpret the Einstein equations, 
written in the form T^ t + T g f av = 0, as a 'micro-balancing' of local stress-energy-momentum 
at each spacetime point: wherever a component of matter stress-energy- momentum is posi- 
tive, the corresponding component of gravitational stress-energy-momentum is negative, and 
vice versa, such that the total is always zero. The idea that —G ab / (8tt) is the local energy- 
momentum tensor of gravity is, of course, a very old idea, first put forward by Lorentz and 
Levi-Civita. It was rejected by Einstein, since it implies that the total energy of a closed 
system would always be zero, which is obviously problematic.^ It is only with hindsight that 
we now realize why the problem was not resolved much sooner. People then did not think 
about boundary terms as much as they do today. Thanks to Brown and York we now know 
that what comes to the rescue is the boundary term in Eq. (|1.5| ). 

In this equation the spacetime is assumed to be the topological product of a three-space 
£ and a real line interval. The boundary component B is a timelike 'tube', topologically 
the product of S — <9£ and the real line interval. (The two spacelike 'end-cap' bound- 
ary components of DM have been omitted, as they play no role in this discussion.) The 
quantity —^—7 IP 6 / (167r), constructed in the usual way out of the extrinsic curvature of 
B, is the gravitational momentum conjugate to the three-metric 7^ induced on B. Now, 
in the spirit of identifying the energy-momentum tensor as the functional derivative of the 
action with respect to the metric, one reads off from Eq. ( |1.5| ) the energy-momentum tensor 
Tg b := — n ab /(87r), which is inherently associated with the boundary B, rather than the bulk 
spacetime. Like any acceptable energy-momentum tensor, TJS b is covariantly conserved (as 
a tensor residing in B) — this follows from the analogue of the diffeomorphism constraint of 
general relativity for the three- surf ace B (rather than E). Physically, B is to be thought of 
as the congruence of world lines of a two-parameter family of observers with four-velocity 
u a , hypersurface orthogonal to a one-parameter foliation of B by spacelike two-surfaces with 
the topology of S. Within their three-dimensional spacetime (B, 7) , the observers measure a 
spatial energy density T^ b u a Ub, which is precisely —k/(8Tr), and thus one is led to Eq. (|1.1|). 
Finally, observe that one can add to the action any covariant functional of the boundary 
three-metric 7^ without affecting the previous argument. This is the source of the reference 
point ambiguity CQE rcf in Eq. (|LJ. This summarizes the central idea of the Brown- York 
approach 0. 

Now if energy is really quasilocal, and calculated via a surface integral involving Tg b , one 
comes to the conclusion that a priori neither T® b . At nor Tg r b av has anything to do with energy! 
While this might be unsettling at first, it is reassuring to know that a satisfactory notion of 
local matter energy density can be recovered from the small sphere limit of quasilocal energy. 
For example, in Ref. [|nj] it is shown that, for a certain choice of reference term CQE ref , the 
Brown- York quasilocal energy contained in an infinitesimal sphere of proper radius r is the 
volume of the sphere (47rr 3 /3) times the local matter energy density T^ t u a Ub (evaluated 
at the center of the sphere) that would be measured by an observer with four-velocity u a . 
Moreover, this is a well established property of most quasilocal energy definitions |19|,|20|, 



4 See the historical discussion given on pages 176-7 in Ref. fTI} ). I thank L. de Menezes for bringing this 
reference to my attention. 
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21, 22, 23L|23], so the result is quite robust. And at higher order in r, gravitational energy 
begins to appear, as will be discussed in detail later. The point is there is no contradiction 
between (i) the local energy-momentum tensor T^ t + T^ av being zero, and (ii) there being 
nonzero stress-energy-momentum in a finite spatial volume. This is because T^ t + T^ av 
is not a local energy- momentum tensor — indeed, if we accept the previous argument, there 
is no such thing. There is only Tg b , associated with the fact that energy is fundamentally 
quasilocal. 

The main purpose of this Introduction is to emphasize, firstly, that energy is fundamen- 
tally quasilocal, i.e., associated with closed spacelike two-surfaces — not points — in spacetime; 
and secondly, there are strong reasons to believe that —k/(8ir) is the correct quasilocal en- 
ergy surface density. The major unresolved problem is how to choose the right well defined 
energy reference term, CQE ref . Rather than address this problem per se, I will begin with 
—k/(8Tr) as an energy surface density and construct a new definition of quasilocal energy 
based on analogy with the special relativity formula: E 2 — p 2 = m 2 . The new definition 
is both physically and geometrically natural, and lies somewhere between the Brown- York 
CQE and the Hawking [25| or Hayward || definitions. A reference subtraction procedure 
is still required, that involves a reference embedding, but this is a codimension-two embed- 
ding that is not subject to the problem that afflicts the Brown- York embedding prescription. 
Moreover, there is a shift in the physics: the reference embedding is not associated with 
determining a reference energy so much as a 'reference angular momentum', so to speak. 
Why angular momentum? Because angular momentum contributes to energy, and the new 
definition can be seen as a precise formulation of this fact at the quasilocal level. 

The paper is organized as follows. In Sec. 2 we introduce the geometrical quantities 
we will use later. Sec. 3 contains the physical and geometrical motivations behind the new 
definition of quasilocal energy (as well as the definition itself). The reference subtraction 
term is discussed in Sec. 4. In Sec. 5 both the spatial and future null infinity limits of the 
energy are examined; the small sphere limit is considered in Sec. 6. Finally, in Sec. 7 we 
examine the new energy in the context of asymptotically anti-de Sitter spacetimes. Since this 
is a lengthy paper I have provided a summary of its results at the end, which also includes 
some additional discussion. 



2 The geometry of two-dimensional spacelike submanifolds 

Let (M, g) be a four-dimensional Lorentzian geometry with signature +2, and S be a 
closed two-dimensional spacelike submanifold. Let u a and n a be timelike and spacelike unit 
normals to S that are orthogonal to each other: u a u a = —1, n a n a = 1, and u a n a = 0. We 
will assume that S is orientable, and an open neighborhood of S in M is space and time 



orientable, so that u a and n a are globally well defined ||26|| . u a and n a are fixed up to an 
arbitrary local boost transformation: 

u' a = u a cosh A + n a sinh A 

n' a = u a sinh A + n a cosh A. (2.1) 
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The physical picture to keep in mind is that of a finite spatial volume E, i.e., a three- 
dimensional spacelike submanifold, whose boundary is S. Although S need not be connected, 
nor simply connected, we will often think of S as having the topology of a three-ball, and S 
that of a two-sphere, and thus will sometimes refer to the direction of n a (assumed outward 
directed) as the 'radial' direction. Given such a three-surface S spanning S it is natural to 
choose u a to be orthogonal to £ (and future directed), in which case n a is tangential to S. 
Physically, u a is the instantaneous four-velocity of a two-parameter family of observers on 
S. With u a thus tied to the spanning surface £, a deformation of £ (preserving S) will in 
general effect a radial boost, Eqs. (|2.1|) . 

The remainder of this section is a summary of some standard facts about the geometry 



of the submanifold S, as can be found, e.g., in Ref. [27], except here we follow a notation 



similar to that used in Ref. [g(J. The surface projection operator, V b , is a tensor defined on 
S by 

V a h := 5 a b + u a u b - n a n b . (2.2) 

(All raising and lowering of indices a,b,c, . . . will be effected with the metric g a b, or its inverse, 
g ab .) A 'surface tensor' is defined as a tensor on S that is left invariant under projection of 
all its indices with the surface projection operator. Obviously one such tensor is the spatial 
two-metric 

(Tab ■= V c a V b g cd = g a b + u a u b - n a n b (2.3) 
induced on S. Another is the corresponding volume form on S, given by 

e«fe := e abcd u c n d , (2.4) 

where e a b c d is the volume form on M. The symbol dS will be used in place of e a b as the 
integration measure for (S,a). 

If V a denotes the Levi-Civita connection of (M,g), then T> a , the Levi-Civita connection 
induced on (S, a), is defined by 

V a Tt = V d a V\V{ ■ ■ ■ V d T e fZ , (2.5) 

where T^" is any surface tensor. Then for any two surface vector fields X a and Y a , the 
Gauss formula reads 

X a V a Y c = X a V a Y c + h c ab X a Y\ (2.6) 

where h c ab is the second fundamental form. Its first index is normal to S, i.e., V d h c ab = 0, 
whereas the remaining two are surface tensor indices, that are symmetric under interchange 
(as can be easily seen by interchanging X and Y in Eq. ( |2.6|) and subtracting the two 
equations). Thus the second fundamental form can be decomposed into components along 
the two unit normals: 

h c ab = u c l ab ~ n c k ab , (2.7) 
where the two extrinsic curvatures are (symmetric) surface tensors given by 

iab = -u c h c ab = v:ry cUd 

k ab = -n c h c ab = V c a VtV c n d . (2.8) 
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It is useful to decompose the extrinsic curvatures into trace and trace-free parts: 

lab = T^Vab + 

kab = ^ka ab + k ab , (2.9) 

where I = <J ab l ab and k = <r ab k ab , and a tilde appearing over any quantity in this paper will 
always mean 'trace- free part of. The mean curvature vector is then 

H c := l -a ab h\ b = ^(lu c -kn c ), (2.10) 

and H ■ H = (k 2 — Z 2 )/4 is the square of the mean curvature. 

Extrinsic curvature is a measure of how a unit normal vector rotates as it is parallelly 
propagated tangent to S in the ambient space (M, g). Two normal vectors means two extrin- 
sic curvatures. However, from Eqs. ( |2.8| ) we see that l a b and k a b measure only the components 
of this rotation tangent to S. There is also a normal component, i.e., the component of the 
rotation of one normal vector along the other. Thus a complete characterization of the 
extrinsic geometry of S requires also the surface one-form 

A a := V b a n c V b u c . (2.11) 

This is an 5*0(1, 1) connection in the normal bundle of S, and its associated curvature 
two-form is 

T ab := V a A b - V b A a . (2.12) 

We will see later that the curvature of the normal bundle of S plays a key role with regard 
to angular momentum. 

While the second fundamental form (including H c and H ■ H) and the curvature of 
the normal bundle are invariant under local radial boosts, the extrinsic curvatures and the 
connection on the normal bundle are not. They transform as 

I'ab = l ab cosh A + k ab sinh A 
k' ab = l ab sinh A + k ab cosh A 

A' a = A a + V a \. (2.13) 

Observe that A a is different from the other extrinsic curvatures in that its transformation 
law, which is a gauge transformation, depends on the derivative of A. 

Our sign conventions are such that the Riemann tensor of (M, g) is defined by (V a Vf, — 
V b V a )X c = R abc d X d} and similarly that of (S,a) by (V a V b - V b V a )X c = K abc d X d (X c is 
a surface one-form in the latter case). Appropriate projections of the Riemann tensor of 
(M, g) yield the Gauss equation: 

VlV{VlV h d R eSgh = TZ abcd + (l ac l bd - l bc l ad ) - ( k ac k bd - k bc k ad ) , (2. 14) 
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the Codazzi equations: 



VlVtV 9 c u h R eSgh = (V a l bc - V b l ac ) - (A a k bc - A b k ac ), 
VlVlV g c n h R eSgh = (V a k bc - V b k ac ) - (AJ bc - A b l ac ), 



(2.15) 



and the Ricci equation: 



V e a Viu 9 n h R, 



(2.16) 



These are the integrability conditions for the isometric embedding of (5*, cr) into (M, g), and 
so by definition of S are necessarily satisfied. 

3 The invariant quasilocal energy 

A physical interpretation of the various geometrical quantities introduced in the previous 
section can be given as follows. The expansion k measures the fractional expansion of the 
area of a small element of S when each point in the element is projected a unit distance 
radially outward. It will have a certain positive value if, for example, S is a round sphere 
enclosing a volume of flat 1R 3 . (For our present purposes, imagine 'flat R 3 ' as a t = constant 
surface in Minkowski space.) Now if 5* is a round sphere of the same area enclosing some 
matter, then, according to the Einstein equations, the matter curves the space inside S in 
such a way that its volume is greater than one would infer by measuring just the area of the 
sphere and using Euclidean geometry. Thus the expansion measured at 5* must be smaller, 
i.e., the areas of spherical shells at larger radii will not increase as rapidly as expected. 
So we see that the unreferenced Brown- York quasilocal energy in Eq. ( |1.1| ) is greater (less 
negative) when S contains matter, than when it does not. This is an intuitive reason why 
A; is a measure of the energy inside S. (It also explains the need to subtract off a reference 
energy of the form given in Eq. ( |1.2|) : k is the nonzero value of k when S merely encloses 
a volume of flat 1R 3 , i.e., no energy.) 

/ is similar to k, except that it measures the expansion of S in time, i.e., in the direction 
of the observer's four- velocity u a . Intuitively, if the observers tend to be moving radially 
outward then the area of the two-surface they are on will be expanding, i.e., / > 0. Con- 
versely, a radially inward motion corresponds to / < 0. Thus I (more precisely, 1/(8it)) can 
be interpreted as a 'radial momentum surface density' p8[ . In the case that (k 2 — I 2 ) is 
positive, the observers can always make appropriate local radial boosts such that I — at 
each point of S, a situation corresponding to a 'quasilocal rest frame'. I will comment on 
this notion more precisely at the end of this section. 

The trace-free quantities k ab and l ab measure the shear of S, and are intimately connected 
with angular momentum (or at least l ab is). For example, consider a set of locally nonrotating 
observers who at coordinate time t are on a constant r, t sphere of the Kerr black hole in 
Boyer-Lindquist coordinates. Their four-velocity is given by 




(3.1) 
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where N is the lapse function, and u)(r,8) = —gt^/g^ is an observer's angular velocity 
as measured from infinity J29|. Starting at Eqs. (|2.8| ) it is not difficult to show that in 
this case / = 0, so here is an example of observers in a 'quasilocal rest frame' as defined 
above. Furthermore, one can show that the nonvanishing components of the shear in the 
time direction are given by 

l H - he - m - B§ . (3.2) 

Physically, a nonzero u reflects the frame dragging caused by the rotating black hole. The 
fact that the degree of frame dragging depends on 9 is what makes the observers at different 
latitudes of the sphere rotate at different rates relative to 'the distant stars', and more to 
the point, relative to each other. This causes a shear effect between observers at neighboring 
latitudes, which obviously disappears when the angular momentum is zero. 

Furthermore, let the locally nonrotating observers label themselves with coordinates 
(#',<//), which at t — coincide with the Boyer-Lindquist (6,(j)) coordinates on 5*. Then 
although the observers always measure the same two-geometry of S as time t goes on, the 
components of the two-metric a a b in their (6' , 0') coordinates will differ from those in the 
(6, 4>) coordinates by a t- and ^-dependent diffeomorphism along the ^-direction. So al- 
though one usually associates shear with a geometrical deformation, for instance a round 
sphere evolving into an ellipsoid, one can also have a physically meaningful shear associated 
with a continuous parameter family of isometric surfaces. This fact plays an important role 
in understanding certain embedding equations we will encounter later, and will be discussed 



in detail elsewhere [30 



The last geometrical quantity to interpret is A a , the connection in the normal bundle. In 
the Brown- York analysis the quantity — A a /(8ir) is called the 'momentum surface density', 
and is denoted as j a 0. The momentum vector j a is tangential to 5*, corresponding to 
a 'rotating two-surface', and thus should be associated with angular momentum. Indeed, 
this is correct: Let B denote the timelike three-surface that is the congruence of world lines 
belonging to the two-parameter family of observers on a two-sphere S. If B admits a Killing 
vector field a , whose orbits lie in S, then one can define the angular momentum charge 

J:= [ dS<p a Ja , (3.3) 
•J s 

which can be shown to coincide with the ADM angular momentum at infinity for asymptot- 
ically flat spacetimes [[?]. Thus we expect both the shear and the connection in the normal 
bundle to play a role in angular momentum at the quasilocal level, and indeed we will see 
that this turns out to be the case. 

Now the first goal of this paper is to provide a physical motivation for the general rela- 
tivists analogue of the special relativity formula: E 2 — p 2 = m 2 . First of all, this formula 
applies strictly to point particles (as opposed to extended objects). One imagines determin- 
ing, say, the instantaneous three- velocity of such a particle by measuring its location in space 
at two closely separated points in time, in some inertial reference frame. In the spirit of the 
quasilocal idea, the analogue of this in general relativity would be to first replace measure- 
ments at a point with measurements on a closed spacelike two-surface S. But measurements 
of what? It would seem that 'measurements of the location of the point particle' 'in some 
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inertial frame' is to be replaced with 'measurements of the two-geometry of S" 'in a generic 
spacetime'. These measurements are to be repeated at two closely separated points in time. 
In the point particle case this yields the three- velocity (or the three-momentum p if one also 
knows m); in the two-surface case it yields l a b, the time component of the extrinsic curvature 
of S. Now I pointed out above that the trace of l a b — more precisely l/(8n) — indeed has the 
interpretation of a momentum: it is the normal (or radial) momentum surface density p8l . 
So it seems reasonable to replace p with l/(8n). What about the trace-free part of 1^1 It 
was argued above that l a b is associated with angular momentum. Insofar as angular mo- 
mentum is qualitatively distinct from linear momentum, its role at least at this point of the 
argument is not clear, and we will simply drop it for now (however, its role will become clear 
later). Notice that dropping l a b is at least roughly consistent with being interested only in 
the two-geometry of S, i.e., the two-metric a a b modulo diffeomorphisms, since, as indicated 
above, l a b is in some cases associated with just diffeomorphisms of S. 

Thus, in the expression E 2 — p 2 we will replace p with 1/(8tt). What should replace El 
Given the previous discussion, the obvious answer is the Brown- York energy surface density, 
—k/(8ir). Clearly 1/{8tt) and —k/(8ir) are on exactly the same geometrical footing, being 
proportional to the timelike and spacelike components of the mean curvature vector H c (see 
Eq. ( |2.10|) ). Thus we arrive at the generalization 

E 2 -p 2 ^ 7^y 2 (k 2 -l 2 ). (3.4) 

Now before we accept this generalization, let us observe that there is something unex- 
pected about it. The four-momentum (E,p ) has become a two-momentum, (— k, 1)/(8it). 
What happened to the other two components of spatial momentum? 1/{8tt) is just the ra- 
dial component; shouldn't the Brown- York momentum surface density j a (in our notation, 
— A a /(8ir)), which is tangent to S, be the analogue of the two missing components of p ? If 
so, then instead of Eq. (|3.4j ) we should have 

E 2 - p 2 ^J^y( k2 ~ l2 - AaA *)- (3-5) 

At first sight this expression is appealing because it manifestly includes a contribution 
from angular momentum, and it is known that in general relativity angular momentum 
contributes to mass. A simple example that illustrates this phenomenon is the Kerr black 
hole, where the ADM mass in excess of the irreducible mass is due to rotational energy. The 
precise relationship is [PJ 

M * = M2 -(50< (3 ' 6) 

where M is the ADM mass, M; r the irreducible mass, and J the angular momentum of 
the black hole. Comparing the right hand sides of the previous two equations suggests we 
conceptually identify |A|/(87r) with the angular momentum term, J/(2M ir ), which seems 
reasonable. This leaves \/k 2 — I 2 / (8ir) to be interpreted as an object like M, viz., a total 
mass, 'total' in the sense that it includes the contribution from angular momentum. But here 
then is the point: \/k 2 — I 2 / (8tt) somehow implicitly already includes the angular momentum 
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contribution to mass. Precisely how will become clear later, but to see immediately that 
this is at least plausible, consider the case I = 0. Then —\/k 2 — I 2 / (8-7r) reduces to the 
Brown- York energy surface density, at least when k is nonnegative, and it is known that the 
(referenced) Brown- York quasilocal energy yields the ADM mass at spatial infinity, which 
includes the correct angular momentum contribution to mass. So we do not need the A a A a 
term in Eq. (|3.5|). Besides, putting it in is counter to our goal of seeing if general relativity 
admits an analogue of the invariant mass, m: While the combination k 2 — I 2 is invariant 
under radial boosts,^ A a A a is not — see Eqs. ( |2.13| ). So from this point of view the right hand 
side of Eq. (|3.5|) is defective, not to mention the generally unsavory fact that it mixes objects 
with different transformation properties.^ The question of 'missing momentum components' 
can also be thought about as follows. A point particle has three components of spatial 
momentum. Likewise, each point on a two-surface S also has three components of spatial 
momentum (more properly, momentum surface density): one normal, and two tangential 
to S. But being tangential, the latter two are associated with a 'rotating surface', and 
hence with angular momentum. In going from a point to a two-surface, two components 
of the linear momentum have become angular momenta. So they do not (directly at least) 
contribute to the expression for an invariant mass given on the right hand side of Eq. ( |3.4| ) 
because, as claimed, this expression already inherently includes the contribution from angular 
momentum. 

Thus we are led to propose the following definition of an 'invariant quasilocal energy' (or 
IQE): 

IQE = -— / alS Vk 2 -I 2 - IQE ref , (3.7) 
8n Js 

where IQE ref is a reference subtraction term that will be defined later. The word 'invariant' 
in 'IQE' refers to the fact just mentioned, that k 2 — I 2 is invariant under local radial boosts 
of the observers on S. And the word 'energy' is used instead of 'mass' — despite our analogy 
between \/k 2 — I 2 / (8tt) and the mass m — because, as we will see in Sec. 7, the IQE behaves 
more like an energy than a mass. So the IQE can be thought of as the amount of 'rest 
energy' contained in S, a quantity independent of the motion of the observers measuring it. 
Notice that the unreferenced IQE is negative. Nominally the reference energy IQE ref is more 
negative, so that the referenced IQE is positive. 

It is useful to express the integrand of the unreferenced IQE in two other equivalent 
forms. Define the pair of null normals £± := u a ± n a on S, and the corresponding null 

5 This was first noted in Ref. [^8). A further discussion of boosted observers in the Brown- York framework 
appears in Ref. [ ^T| . 

6 Hayward's Q definition of quasilocal energy includes an angular momentum contribution of the form 
— Lu a uj a , analogous to the —A a A a term in Eq. ( |3.5| ). Hayward's uj a is a suitably normalized anholonomicity, 
or 'twist', of the pair of null normals to S, and encodes essentially the same information as A a . The important 
distinction is that, unlike the connection A a , the object u> a is boost invariant, and so representing angular 
momentum with a term proportional to —Lo a uj ai as Hayward does, is perfectly acceptable. (The relationship 
between A a and uj a is discussed in Appendix B of Rcf. p^|.) However, there is no need, or even natural 
way for uj a to enter our work here. For instance, the A a A a term in Eq. (3.5) cannot simply be replaced 
with uj a uj a , since the (tentative) inclusion of this term is suggested by the physical interpretation of A a as a 
momentum surface density. This interpretation arises from Brown and York's Hamilton- Jacobi analysis of 
the gravitational action 0, and it is not clear that a similar interpretation can be given to u> a . 
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expansions 



9±:=o*V a £±b = l±k, (3.8) 
cf. Eqs. fl2.8| ) and ( |2.9| ). Then we have the following three equivalent expressions: 



1 — _ l r— 1 



— Vk^P = ■jr-J-O+O- = j-VWTh. (3.9) 

87T 07T v 47T 



For the last expression recall the definition of the mean curvature given after Eq. ( 2.10|) . 



Thus the unreferenced IQE in Eq. (|3.7| ) has a very simple geometrical interpretation: up to 
a proportionality constant, it is just the mean curvature of S, averaged over S. Because of 
the square root it is defined only when k 2 — I 2 > 0, i.e., at each point of S the mean curvature 
vector H c in Eq. ( |2.1U| ) must be either spacelike or null, but never timelike. Roughly speaking, 
this means that the area of S changes more rapidly in a radial direction, than in time. For 
example, this condition is satisfied for the constant r, t two-spheres outside the horizon of a 
Schwarzschild black hole, but not for those inside; on the horizon the unreferenced IQE is 
zero. 

In terms of the null expansions, recall that a future (past) trapped surface is one for which 
both ingoing and outgoing null expansions, 6L and 6 + , are everywhere negative (positive) on 
S KJ. Thus, the unreferenced IQE is imaginary when S is a future or past trapped surface. 
It is real only when no point on S is 'trapped'. Now a future trapped surface does not quite 
characterize a black hole, and more subtle characterizations have been proposed for a local 
definition of a black hole horizon For example, Hayward [^] has introduced the 

notion of a 'future outer trapping horizon', H, characterized by: (i) #_|# < (in-going light 
rays converging), (ii) 9 + \h = (outgoing light rays instantaneously parallel on the horizon), 
(hi) 9 + \h+ > (outgoing light rays diverging just outside the horizon), and (iv) 9 + \ H - < 
(outgoing light rays converging just inside the horizon). According to this general definition 
of a black hole, the unreferenced IQE is nonzero just outside the horizon, zero on the horizon, 
and undefined (or imaginary) just inside the horizon. In this connection see also Ref. [ |33|j . 

Furthermore, observe that the condition for the integrand of the unreferenced IQE to be 
real and nonzero, namely k 2 — I 2 > 0, is precisely the same condition that ensures that the 
observers can always, by appropriate local boosts, go to a 'quasilocal rest frame' in which 
I = at each point of S. Such a two-surface is analogous to a massive particle. The case 
k 2 — I 2 = everywhere on S, for instance when S is a future outer trapping horizon, is 
analogous to a massless particle, for which no quasilocal rest frame exists. And finally, the 
case k 2 — I 2 < 0, say 'inside' a future outer trapping horizon, corresponds to a superluminal 
particle.[] 

The situation is actually more subtle than indicated in the previous two paragraphs. The 
conditions for ^—9 + 9- to be real are reminiscent of the condition < required for the 
holomorphic case of the Dougan-Mason quasilocal energy to be nonnegative. The conditions 



7 As noted in the text, the spirit of the quasilocal idea is to replace measurements at a point (of certain 
aspects of a point particle, say) with measurements on a closed spacelike two-surface. If one takes seriously 



that Eq. (3.4) is the generalization of point particle rest mass, then one is quickly led to speculate that a 
closed spacelike two-surface is the generalization of a point particle. This is curiously reminiscent of string 
theory, except that the one-dimensional string is replaced by a two-dimensional surface. 
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9 + > and 6L < essentially imply that the two-surface S is 'suitably' convex |TjJ. To 
emphasize that " 'suitably' convex" is not a serious restriction, in particular it does not mean 
that S cannot be concave, consider a two-parameter family of observers at rest in an inertial 
frame in flat spacetime. Suppose that at t — they lie on a two-sphere S that is round 
except for a small indentation. Then I = at each point of S, and k is positive everywhere 
except in a small region near the center of the indentation, where it is negative. Thus there 
will be a circle of points C at which k = 0. So at each point of S we have k 2 — l 2 > 0, equality 
holding on C. One might worry that a radial boost at a point on C will make I 2 > 0, and 
hence k 2 — I 2 < 0. But of course this will not happen: If we consider a second set of observers, 
boosted relative to the first, then k = I = on C implies k' — I' — on C. So we can consider 
the second set of observers to be boosted radially outward in the region of the indentation, 
such that the indentation, and its attendant set of fixed points C, smoothly disappear as the 
sphere evolves in time, k switches from negative to positive by passing through the origin 
of a k-l diagram. Thus we can imagine a wide class of two-surfaces, including ones with 
indentations, and dynamically changing in time, for which k 2 — I 2 > everywhere on S. 
Moreover, bear in mind that the observers are allowed to accelerate, so there is a great deal 
of freedom for them to maintain a 'physically reasonable' S. Nevertheless, what is needed 
here is a careful analysis based on Raychaudhuri-like equations for a two-parameter family of 
accelerated timelike curves (as opposed to the more usual case of timelike geodesies). Such 
a detailed analysis is outside the scope set for this introductory paper. 



4 The reference invariant quasilocal energy 

As in the Brown- York case, the unreferenced IQE diverges in an asymptotically flat 
spacetime as the two-surface S is taken to (spatial or null) infinity, and so must be regulated 
with a reference term, IQE ref , as already anticipated in Eq. (|3.7|). To better understand 
the nature of our definition of the invariant quasilocal energy, and to help suggest a natural 
choice for IQE ref , we now make use of the Gauss embedding equation given in Eq. ( |2.14j ). 



This equation has only one independent component. Transvecting both sides with o ac o M 
reduces it to the scalar equation 

ooR = 1l- X -{k 2 -l 2 ) + (~k 2 -P), (4-1) 

where aaR is shorthand for a ac a bd R a b c d, {k 2 — I 2 ) is shorthand for (k ab k a b — l ab lab), and 1Z is 
the scalar curvature of (S, a). Using this equation we can express the IQE given in Eq. ( ^.7| ) 
in the equivalent form 

IQE = -— / dS J 2 111- aaR + (k 2 - I 2 )] - IQE ref . (4.2) 
8n Js V L J 

We remark here that aaR is a natural geometrical object called the sectional curvature of 
(S, a) as embedded in (M,g) It will play an important role in what follows. 

Now the definition of the unreferenced IQE is rooted in the extrinsic geometry of the sub- 
manifold (S, a), thought of as a two-surface isometrically embedded in the spacetime (M, g). 
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It is then natural to define the reference IQE to be of the same form as the unreferenced IQE 
in Eq. Q4.2| ), i.e., to be the same geometrical object, except with (S, a) now isometrically em- 
bedded in a different spacetime — some reference spacetime, (M ref , g rei ). Thus IQE ref will be 
the integral in Eq. ( |4.2j ) (or ( |3.7| )), except with all quantities referred to the reference space- 
time, which we indicate with a superscript 'ref. Note that although the extrinsic geometry 
of S will be different in (M ref , g rei ), its intrinsic geometry, by assumption, will not. So in the 
IQE ref integral we are constructing we can set dS r ^ = dS, lZ Tei = 1Z, and (aaR) rei = aaR mi . 
Also note that, in general, M ref ^ M (topologically). For example, S may be a two-sphere 
embedded in a black hole spacetime, with M = R 2 x S 2 , whereas the reference spacetime 
might be Minkowski space, with M ref = R 4 . With this understanding, we define 

IQ E ref = --L f dS J(k2-pyef = -l- [ dS I2 \n - aoBT* + (k 2 ~ l 2 Y e{ ] ■ (4.3) 

07T JS o7T JS V L J 

The term aaR rc{ is shorthand for (y ac a bd R r ^ cd) where R a e b { cd * s the Riemann tensor of the 
reference spacetime. 

Typically one is motivated to choose a reference spacetime of constant curvature, the geo- 
metrical reason being that then the Gauss, Codazzi, and Ricci embedding equations make no 
reference to 'where' (S,a) is embedded in (M ref , g ref ). In other words, the conditions placed 
on k a f, l a f, and A a ef by the 'reference version' of Eqs. ( |2. 14^pTT6|) — which are just integra- 



bility conditions for the reference embedding — do not depend on knowing the embedding 
itself |[34| . This is a pleasing criterion because it keeps the reference spacetime 'abstract', 
rather than 'concrete'. For a four- dimensional space of constant curvature we have 

Robed = -^(fahbd - 9bc f 9ad)i (4-4) 

where C is the constant value of its scalar curvature. For this choice of reference spacetime 
one gets 

aaR rci = -. (4.5) 
6 

For example, for Minkowski space we have C = 0, and for anti-de Sitter space we have 
C = —12/£ 2 , where £ is the radius of curvature of the anti-de Sitter space, and is related to 
the (negative) cosmological constant A by A = —3/£ 2 . We will return to these two examples 
later. 

The idea of embedding (S, a) into some reference space(time) is in the same spirit as 
the Brown- York approach, but an important difference that arises out of using the invariant 
quantity \/k 2 — I 2 , rather than k, deserves further comment. In the Brown- York approach k 
is the trace of the extrinsic curvature of (S, a) as embedded in a three-geometry, (E, h), where 
S is a spacelike three-surface spanning S, with induced metric h a ^. Thus it is natural to take 
k mf as the trace of the extrinsic curvature of (S, a) as embedded in some three-dimensional 
reference space, (S ref , h Tei ). So the embeddings of 5*, for both the unreferenced and reference 
CQE, inherently have a t/iree-dimensional target space. On the other hand, \Jk 2 — I 2 is 
proportional to a geometrical invariant of S, namely its mean curvature, and makes no 
essential reference to a spanning three-surface E (making the IQE 'truly quasilocal' in the 
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sense that it depends on S aloneQ). As a consequence, the embeddings inherently have a 
/owr-dimensional target space(time). 

The advantage of a three-dimensional target reference space, say flat R 3 , is that when 
the embedding exists, it is unique (up to translations and rotations), and so the Brown- York 
CQE ref is unique. The disadvantage, as is well known, is that such embeddings do not exist 
for all (S, a) of interest, and this problem is not limited to just a few isolated exceptional 
cases. 

For a four- dimensional target reference spacetime, say Minkowski space, the situation 
is reversed: an embedding always exists, but it is not unique. Regarding the first half 
of this statement, Brinkmann |J5| has shown, by a simple explicit construction, that any 
n-dimensional conformally flat Riemann space can be considered as a particular cut of a 
light cone in (n + 2)-dimensional Minkowski space. And conversely, any cut of such a light 
cone gives an n-dimensional conformally flat Riemann space. Now any n = 2 space is, of 
course, conformally flat, and thus any (S, a) can always be so embedded, even if S has 
regions with negative scalar curvature. In the Introduction I mentioned the example of 
the horizon of the Kerr black hole, which cannot be globally embedded into flat R 3 when 
the angular momentum exceeds the irreducible mass, which coincides with the two-sphere 
developing regions with negative scalar curvature |16 |. However, it is a simple exercise to 



apply Brinkmann's construction and thus globally embed the horizon into a light cone of 
four-dimensional Minkowski space. I will omit the details of this calculation, and just note 
that the embedding is valid for all angular momentum J (up to and including the extremal 
black hole case), and changes smoothly with J, including at the critical point when J equals 
the irreducible mass. 

On the other hand, in a co dimension-two (versus -one) embedding there is more 'elbow 
room', and consequently the embedding is not unique — there is a function worth of freedom 
(which will be discussed in detail in Ref. 0]). This results in an ambiguity in the reference 
energy, IQE ref , which enters via the reference shear term {k 2 — Z 2 ) ref — see Eq. (|4.3| ). This is 
the only term in IQE re not yet determined, and the only one for which we require a reference 
embedding. Observe that in the Brown- York approach the undetermined quantity is k rei , an 
expansion. Here it is (k 2 — Z 2 ) ref , a shear. I will now argue that it is precisely this term that 
plays the key role in properly incorporating angular momentum into the IQE. The basic idea 
is simple, but first we will introduce some notation. 

In the previous section we introduced the null normals £± = u a ± n a , and corresponding 
null expansions 9± = I ± k in Eq. fl3.8|) . Similarly, the (trace-free) shears in the two null 
directions are defined by s± a b := l a b ± k^. The curvature of the normal bundle has only 
one independent component, and can be written as T a b = (3 r /2)e a b for some scalar field 

where e a b is the volume form on S defined earlier in Eq. ( |2.4j ). With this notation, and 
assuming that the reference spacetime is one of constant curvature, i.e., Eq. ( f4.5|) holds, the 

8 The CQE can also be made 'truly quasilocal', in a slightly different sense: by relaxing the restriction that 
the foliation of the spacetime (i.e., E) be orthogonal to the boundary B M, it is shown in Ref. J3l| that the 
resulting CQE no longer depends on E, but instead just depends on the foliation of B. 
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Gauss, Codazzi, and Ricci embedding equations given at the end of Sec. 2 take the form 



j = n + \efe^ - s f v_ f b a , (4.6) 

= \(V a t Af)df - (V b =f AfKi 6 a , (4.7) 
= + \e\[sf , s^] b a . (4.8) 

In the Ricci equation, [s+ f ,s^ f ] 6 a denotes the commutator of the shears: s r f b c s r ^ ic a — 
s ref fe cS ^e f c a _ Notice that by using null directions, rather than u a and n a , the Codazzi equations 
have decoupled into a '+' and a ' — ' set. 

Our task is thus: Given a ab , and hence TZ, V a , and e ab , solve these embedding equations 
for the unknown quantities 8± f , s± f ab , and A^- (Of course J-" ref = 2e ab T) a A r b ci is not an 
independent quantity.) In particular, we are interested in the solution for the boost invariant 
reference shear term 

(p-r 2 r f = - s r + ef v r _ ef6 , (4.9) 

appearing in the Gauss equation, which is to then be substituted into the second integral of 
Eq. ( f4.3|) . Or equivalently, solve for (k 2 — Z 2 ) rcf = —B™ 9_ and substitute the answer into 
the first integral of Eq. ( |4.3[ ). This is how IQE ref is determined. 

However, as already noted, any solution we obtain is not unique. We can see this im- 
mediately by counting functional degrees of freedom. 9± are two functions, s T ± ab are four 
(the two shears are symmetric and trace-free), and A r ^ { are two. These eight functions are 
subject to six equations: Gauss is one, Codazzi are four, and Ricci is one. This leaves two 
arbitrary functions in the solution. But owing to the invariance of the embedding equations 
under a local boost transformation (see Eqs. (|2.13|) ), one of these functions is just the boost 
parameter, A, leaving one nontrivial arbitrary function in the solution. 

The question then arises, Is there a natural way to impose one additional functional 
condition on the unknowns so that the embedding, subject to this additional condition, is 
unique, and hence IQE ref is unique? One of the central ideas in this paper is to impose the 
additional condition 

F ci = T, (4.10) 

i.e., the curvature of the normal bundle of 5 as embedded in the reference spacetime 
(M , g rei ) should equal that of 5 as embedded in the original physical spacetime, (M,g). 

There are several reasons why it is geometrically natural to demand T vei = T . First, the 
two-surface S has two connections: one is an 50(2) connection on the tangent bundle of 5, 
associated with the curvature TZ a bcd (which has only one independent component, namely 
TV), and the other is an 50(1, 1) connection on the normal bundle of 5, associated with the 
curvature jF afc (which also has only one independent component, namely J 7 ). In fact both 
of these connections are metric connections, associated with the metrics in the tangent and 
normal bundles to 5, respectively [p7| . Furthermore, Szabados [|(| has considered the two- 



dimensional version of the Sen connection for spinors and tensors on a submanifold such as 
5, and has found that the two-surface spinor curvature has, essentially, imaginary part equal 
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to TZ, and real part equal to T. Finally, although A a is a measure of extrinsic geometry, 
as pointed out earlier it is not really on the same footing as the extrinsic curvatures k a b 
and lab, since its transformation law under local radial boosts is qualitatively different — see 
Eqs. ( |2.13D . It transforms like the connection that it is, and gives rise to a curvature, and so 
arguably has more in common with TZ than with k a b and l a b- The point is, T is really on the 
same geometrical footing as TZ. We have already demanded that TZ rei = TZ, as a necessary 
condition for the embedding of {S,a) into (M ,g rei ) to be isometric. So demanding also 
that F rei = T is thus seen to be quite natural. 

Unfortunately, implementing Eq. (|4.10|) seems like an intractable task. Embedding equa- 
tions involving curvature of the normal bundle, i.e., codimension-two (and higher) embed- 
dings, have, of course, been studied for a long time. With regard to solutions, although one 
expects to be able to express ik 2 — l 2 ) mf in terms of T, TZ, and their derivatives, I am not 
aware of any such general results in the literature. In fact, much of the literature on such 
embeddings considers the case T = 0, which is not the case we are particularly interested in 
here (a notable exception is Ref. [J27|J ). One possible way to proceed is as follows. Given jF rcf 
(= T by Eq. (|4.10|) ), choose A r ^ such that T rei = 2e ab T> a A T b ei . There may be a convenient 
gauge choice, such as 7) ■ A vcf = 0, or Z rcf = 0. Then view the Codazzi equations (O) as a 
set of four linear partial differential equations for the four independent degrees of freedom 
in the two shears. These equations are of second order if one makes use of the fact that 
any trace-free symmetric tensor s a b on a two-surface (S, a) with two-sphere topology can be 
expressed as s a b = 7) a Vb+7)bV a — a a bV-v for some vector field v a . Thus solve for s r ± ab in terms 
of the expansions, 8± f , and their derivatives. The expressions one obtains at this stage are, 
in general, nonlocal. Then substitute these into the Gauss and Ricci equations, (|4~6l ) and 
( |4.8| ), which are really just nonlinear algebraic constraints. But because the shears involve 
nonlocal operators acting on the expansions, one ends up with two nonlocal and nonlinear 
partial differential equations for the two expansions. Remarkably, it is almost possible to 
solve these equations, but in the end one encounters a certain combination of nonlocality 
and nonlinearity that makes the final step to a solution seem impossible. Nevertheless, it 
appears that the solution for (k 2 — / 2 ) ref , if it can be found, almost certainly depends in a 
simple way on both TZ and T , and derivatives (of a finite or possibly infinite order) of these 
two curvatures. I am suggesting that it is through this subtle presence of T in {k — / 2 ) ref 
that angular momentum is properly incorporated into the IQE. 

So although a direct attack on the embedding equations has not yet yielded a solution, 
fortunately one can make some progress of a general nature by calculating the first and second 
order variations of ,F ref and (k 2 — / 2 ) ref under isometric deformations of a given embedding. 
The idea is to see how both of these quantities change under such a deformation, and 
thereby infer how (k 2 — Z 2 ) ref depends on jF ref , and hence angular momentum. The results 
are somewhat involved, and will be given elsewhere JHJ. For now let us start by making 



some simple observations regarding the enigmatic object (k 2 — / 2 ) rcf . 

To begin with, one might object to our argument thus far because it implies that (k 2 — 
l 2 ) ret , and hence IQE ref , depends on the extrinsic geometry of S as embedding in the physical 
spacetime. In particular, through Eq. ( |4. 1Q ) and the reference embedding equations, (k 2 — 
/ 2 ) ref depends on T . On the other hand, it is often stated that a reference subtraction term 
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should be a functional of only the intrinsic geometry of (S, a). However, notice that there 
is no dependence on the extrinsic curvatures proper, i.e. l a b and k a b, only a dependence 
on J 7 , a quantity which I argued above is really on the same geometrical footing as the 
intrinsic quantity TZ. Moreover, as discussed in the Introduction (refer to Eq. ( |1.5| )), in the 
Brown- York approach one is free to add to the action any functional of the boundary three- 
metric, 7 (,, which contains information about the two-metric a a b, as well as information 
about how (S, a) is embedded in the three-boundary B. For instance, one could add to the 
action a boundary integral of the scalar curvature of B, whose variation would add to IP 6 



in Eq. ( p..5|) a term proportional to the Einstein tensor of 7 a fe, as is done in Ref. |36|. Such a 
term obviously depends on some extrinsic geometry of (S, a). In Brown and York's work this 
fact is of course recognized, but being in a Hamiltonian framework, they restrict the form 
of the arbitrary boundary functional such that the energy surface density (— fc/(87r)) and 
momentum surface density (— A a /(87i)) of S in a particular spacelike hypersurface £ depend 
only on the canonical data on E. This effectively means that their reference subtraction 
term can depend only on a ab [0. But as I emphasized earlier, our approach is based on the 
invariant object \Jk 2 — I 2 , and makes no essential reference to a three-surface £ spanning S. 
The invariant quasilocal energy constructed here does not come out of a canonical analysis, 
so there is no reason that our subtraction term cannot depend on T . 

So the shear term {k 2 — l 2 ) rc{ is allowed, but is it really necessary? Perhaps it is just 
an unsavory term resulting from a poor definition of the IQE. For instance, looking at the 
Gauss embedding equation fl4.1|) one might be tempted to write, instead of Eq. ( |3~4l ) . 



E 2 ~ P 2 j^- 2 [(k' 2 - I' 2 ) - 2(k 2 - 1% (4.11) 



where the additional shear term on the right hand side is perhaps the proper way to include 
angular momentum, somewhat like the A a A a term we attempted in Eq. (|3.5| ). This would 
have the advantage of changing IQE ref in Eq. (4.3) to 



IQE rcf = --^ f dS ^2 [R-aaR^], (4.12) 

o7T J S 



which is clearly unique, and moreover, requires no reference embedding. Eq. ( |4 . 1 2| ) is a more 



general case of the zero point energy suggested by Lau |37j (except that his derivation of 



it requires a reference embedding — we will return to this point later). But unfortunately it 
cannot be correct. For example, when the reference spacetime is Minkowski space, Eq. ( |4.5|) 
tells us that aaR vci = and thus the radical reduces to V2TZ, which is not defined for 
negative TZ. Nor is this problem properly solved by taking C ^ in Eq. ([4.5|) , since this 
would put an ad hoc fixed lower bound on TZ. 

So the shear term (k 2 — l 2 ) rei is not only allowed, it is necessary (or at least its absence 
leads to an unsatisfactory result). In fact its role seems to be to keep nonnegative what 
is under the square root in Eq. ( [4.31 ). To see this more clearly, consider the special case 
that (S, cr) is embeddable in flat H 3 . If our reference spacetime is Minkowski space, we can 
then choose to embed (S, cr) in a t = constant slice and, within this slice, the embedding is 
essentially unique. In this case it is easy to see that we will have /Jj| f = 0. And by assumption, 
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aaR vcl = 0, so Eq. QO) reduces to 



IQE ref |, rcf=0 = — / dS\k iel \ = -— I dSJ2 K+(k^) 2 . (4.13) 
ah 8n Js 8n Js V L J 

The uniqueness of the embedding means that k rel and fcjjf are unique. Now clearly, no matter 
what the surface is, the spatial shear must be such that what is under the square root 
in Eq. (}4.13|) is nonnegative, because |fc ref | is real. (More properly, one should look at the 
'reference version' of Eq. ( |4.1| ), with (P ef ) 2 = (/ ref ) 2 = <j<jR tc{ = 0.) For example, consider a 
dumb-bell-shaped surface of revolution in flat R 3 . In a region near the throat of this surface 
1Z is negative, nevertheless at every point of the surface we have 1Z + (k rcf ) 2 > 0. 

I emphasize that, even when (S, a) can be embedded in flat R 3 , its embedding in 
Minkowski space need not be chosen to be in a t = constant slice, as was done in the 
previous paragraph. One may also embed it in a light cone, or in a host of other ways — 
remember that there is a function-worth of freedom in our choice. I argued in the context of 
Eq. (|4.10| ) that this freedom has to do with angular momentum, or more precisely, the cur- 
vature of the normal bundle. For the t = constant embedding, Z^ f = implies s± f afe = ±A£| f , 
and so inspection of Eq. ( |4.8| ) reveals that in this case jF ref = 0. However, it is not hard 
to see that starting with such a t = constant embedding one can perform an infinitesimal 
isometric deformation of the embedding out of the t = constant plane, i.e., in a direction 
ipd/dt, where tp is an arbitrary function. Furthermore, I show in Ref. that after such an 



infinitesimal deformation jF ref is no longer zero, and can in fact be made to be essentially 
any infinitesimal function we like by a suitable choice of (p. It is not hard to imagine (just 
hard to do!) that by integrating such isometric deformations one may be able to achieve a 
two-geometry (S, a) isometrically embedded with any desired curvature of the normal bun- 
dle. With this in mind, it would seem unnatural to reference-embed, e.g., a constant r, t 
two-sphere of the Kerr geometry (which is easily shown to have T ^ 0) as a two-sphere 
in a Minkowski reference spacetime with jF rcf = (say, in a t = constant slice), when it 
seems possible to instead embed it with jF ref = T . Note, however, that as the embedding is 
deformed out of the t = constant surface, will also cease to be zero, and so will introduce 



a negative contribution to the quantity under the square root in Eq. (|4.3|). This jeopardizes 
the nonnegativity of this quantity. But at the same time we clearly cannot simply throw 
away the — (/ ref ) 2 term, since this would violate a key property of the IQE, namely its in- 
variance under local boosts. Short of solving the embedding equations for a generic jF rcf 
and explicitly checking, I do not know of any guarantee of nonnegativity. Equivalently, in 
question here is the nonnegativity of the quantity (k 2 — l 2 ) rei , which is the same type of 
question as the nonnegativity of {k 2 — I 2 ) discussed at the end of Sec. 3. I do not at present 
have a complete answer to either of these difficult questions. 

Finally, one might guess that it is possible to avoid the embedding problem entirely by 
simply setting [k 2 — Z 2 ) ref = (k 2 — I 2 ), which is in the same spirit as Eq. ( |4.1U[ ) in that, like 
J 7 , the shear term has something to do with angular momentum. But this does not work. 
For example, it is not hard to show that, although [k 2 — I 2 ) is in general nonvanishing on 
constant r, t spheres of the Kerr geometry, it happens to vanish on the horizon. So if we set 
(k 2 — Z 2 ) ref = {k 2 — I 2 ), then in calculating IQE ref for the Kerr horizon example we would 
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run into the same problem with negative 1Z as we did in Eq. ( |4.12 ). So such a prescription 
must not be valid, and we cannot avoid the embedding problem this way. 

Let us conclude this section by addressing the question, What is the relationship between 
the IQE defined here and the Brown- York CQE? We begin by supposing that the following 
four conditions are satisfied: (i) (S, a) is a two-surface in the physical spacetime such that 
k 2 — I 2 > 0; (ii) k > 0; (iii) (S, a) is such that it can be embedded in flat 1R 3 ; and (iv) for the 
embedding in (iii), /c ref > 0. As discussed earlier, condition (i) ensures that the unreferenced 
IQE is well defined — roughly speaking, S is strictly outside of a black hole. Then it is always 
possible to go to a 'quasilocal rest frame' where I = on S, and the integrand in Eq. (|3.7|) is 
just |fc|. Given condition (ii), the unreferenced IQE thus reduces to the unreferenced CQE 
in Eq. (|1.1|) , provided the observers in the Brown- York case are in a 'quasilocal rest frame'. 
Condition (iii) ensures that the Brown- York prescription is well-defined, and allows us to 
choose a i = constant embedding in Minkowski space, as above, and get Eq. ( |4 . 1 3| ) . The 
first integral in this equation, together with condition (iv), shows that our IQE ref reduces 
to the Brown- York CQE ref in Eq. flTJ). So if these condit ions hold, and we choose to use 
a t = constant embedding to calculate IQE ref , then our invariant quasilocal energy is the 
same as the Brown- York 'rest energy'. In most applications considered in the literature these 
conditions are satisfied, and the IQE will then share all of the desirable properties of the 
CQE. For example, it will be the thermodynamic energy that appears in the first law of 
black hole thermodynamics for Schwarzschild black holes, as considered in Ref. 0. 

On the other hand, I emphasize that conditions (ii) and (iv) are easily violated. One 
need only think of a round sphere with a small indentation, an example discussed at the 
end of Sec. 3 (except here the spacetime is generic). So in general, the IQE defined here is 
not simply the Brown- York 'rest energy'. Furthermore, we need not choose a t — constant 
embedding to calculate IQE ref . Indeed, as I have argued, such a choice is unnatural when 
T 7^ 0. In short, the invariant quasilocal energy defined here is not quite the same object as 
the Brown- York quasilocal energy. Note that the aforementioned thermodynamic nature of 
the CQE is derived in Ref. |F] assuming that S is a round sphere. It would be interesting to 
extend this analysis to indented spheres, for instance, and determine which, if either of the 
CQE or IQE, is the 'correct' thermodynamic energy. 



5 The large sphere limit of the IQE 

Let us now assume that the spacetime (M,g) is asymptotically flat, and evaluate the 
limit of the invariant quasilocal energy as (S, a) tends to a large sphere at infinity. At 
spatial and null infinity we might expect these limits to be the ADM and Bondi-Sachs 
masses, respectively. Let us see if this is so. 

Let r be a 'time' function on M such that r = defines a spacelike (respectively, null) 
hypersurface TC Tf of topology Rx S 2 extending to spatial (respectively, null) infinity. Letting 
the parameter r* vary over some range gives a foliation of a part of M. Let r be a function 
on M such that r = r* defines a hypersurface that intersects each leaf 7i Tt (over the allowed 
range of r*) in a spacelike two-sphere, S Tt ^ t . The parameter r* ranges to infinity, and over 
its range the surfaces 5V,,r, provide a foliation of Ti Tt . We are interested in the limit — > oo, 
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with r* arbitrary but fixed. In a rather benign abuse of notation we will refer to SV*,r* as 
simply S, and 'take the limit as r — ► oo with r fixed'. The metric induced on S will, as 
usual, be denoted as a a b (in abstract index notation). 

Now assume that the functions r and r have been chosen such that (S, a) tends to a round 
sphere at infinity. Thus the components of its metric in spherical coordinates x l = (9, (ft) 
have an asymptotic expansion of the form 
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Z sin 2 9 



+ 0< 



(5.1) 



In this expansion, X, Y, and Z are each arbitrary functions of r, 9, and <fi. The symbol 
< (r~ n ) denotes a term that falls off faster (or grows more slowly, depending on the sign of 
n) than r~ n , but not necessarily according to a power of r. For example, rather than 0(1), 
the remainder term 0<(r) might grow as lnr. The motivation for this increased generality 
will be explained below when we consider the large sphere limit at null infinity. Furthermore, 
we can choose (the function) r to be an areal radius, in which case we may take \fa = r 2 sin 6, 
where a = det a^. It is easy to see that this requires Z = —X in Eq. ( |5.1|) . 

The scalar curvature of a round sphere of areal radius r is 2/r 2 . Since the metric in 
Eq. ( |5.1[ ) differs from that of a round sphere by a term one power lower in r, we immediately 
have that its scalar curvature 1Z is given by 
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1 + 

^2 /y^> 



(5.2) 



where the remainder term A^, is of order one.0 In our asymptotically flat spacetime the 
components of the Riemann tensor fall off as 1/r 3 , and so the same will be true of aaR, the 
sectional curvature of (S, a). In the present context the appropriate reference spacetime is 
Minkowski space, and so aaR rci = 0. The only other terms to consider in Eqs. ( |4.2| ) and 
( |4~3l) are the shear terms, (k 2 — I 2 ) and {k 2 — / 2 ) ref . We will see below that these, too, fall off 
at least as fast as 1/r 3 in both the spatial and null infinity limits. In the large sphere limit 
the unreferenced IQE thus behaves as 
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(5.3) 



The reference IQE behaves similarly, except we have aaR re{ = 0. Thus 




^-0 + (P-/ 2 ) 
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-2> 



(5.4) 



9 In Ref. j38| it is shown that = T> ■ v + 0<(1) for some vector field v in S. In other words, to leading 
order An is a divergence. This plays a crucial role in some of the results in Ref. J3^]. However, we will 
not need to use this fact, except to provide some insight into our discussion of the 'solution' of the Ricci 
embedding equation in Sec. 5.2. 
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In forming the difference of the previous two expressions it is important to observe that not 
only do the divergent terms coming from the 2/r 2 piece of 1Z cancel, but also the (finite) 
remainder terms are the same in both, and thus also cancel, independent of what A^ 
is. Thus we find that the large sphere behavior of the (referenced) IQE is given by 

IQE = -L / dSr [aoR - (P - I 2 ) + (P - Z~ 2 ) rof + 0( r -% (5.5) 
Io7r Js 

Of course if the shear terms fall off as 1/r 4 or faster they get absorbed into the 0(r~ 4 ) 
remainder term. It is worth emphasizing that in the large sphere limit the square root in 
the IQE is eliminated by the fact that TZ dominates over the other terms. The areal radius 



factor r outside the brackets in Eq. (|5.5|) is really y 2/7^.0 In Sec. 6 we will see a similar 
mechanism at work in the small sphere limit. But in the intermediate regime the IQE is, in 
general, an integral of the difference of two radicals. 

As a quick check of Eq. ( |5.5|) let us evaluate the right hand side for the Schwarzschild 
geometry. In the usual Schwarzschild coordinates r is an areal radius, and it is a simple 
exercise to compute the sectional curvature of a t, r = constant two-sphere. The result is: 
ooR = 4M/r 3 . The shear terms obviously vanish, and with dS = r 2 dQ (dQ the measure on 
the unit round sphere) one immediately gets IQE = M, the ADM mass of the black hole. 



Now the main task is to investigate in detail the shear terms in Eq. (5J5), which we will do 
separately for the spatial and null infinity limits, respectively. 

5.1 The spatial infinity limit 

Rather than proceed with complete generality, it is more instructive to consider an asymptot- 
ically flat metric that exhibits angular momentum explicitly, and then see how this angular 
momentum works its way into the shear terms. (We will be completely general in the more 
interesting null infinity limit case.) The spacetime far from any isolated stationary (nonra- 
diating) rotating source is described asymptotically by the Kerr metric (see Sees. 19.3 and 
33.3 of Ref. f29[), so let us take (M,g) to be the Kerr spacetime. We choose the following 
basis of orthonormal one-forms: 

„0 _ AT A+ „1 



e = N dt, e = —=dr, e = pd6, e = Jg^dcj) — udt), (5.6) 




which are associated with locally nonrotating observers. The corresponding basis of or- 
thonormal vector fields is 



The notation used is standard: x a = (t,r,9,(fi) are Boyer-Lindquist coordinates, u>(r, 




~9t<t>/ 9<jxt> is an observer's angular velocity as measured from infinity, N = Juj 2 g^ — g tt is 



10 This mechanism works even when the sphere is not asymptotically round. In this case the shear terms 
contribute at a higher order, viz. 1/r 2 , in an effort to keep what is under the square root sign positive, as 
discussed earlier. In other words, the factor \/2/TZ is modified in such a way that negative 1Z is likely not a 
problem. We will not consider this more complicated case here. 
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the lapse function, etc. (see, e.g., Sec. 33.4 of Ref. P5|). Let A,B,... be indices labeling 
the basis vectors and one-forms, ranging from to 3, and I, J, . . . denote the subset of these 
taking values 2 and 3. These indices are raised and lowered with the flat Lorentz metric 
Vab = f] AB = diagonal(-l, 1, 1, 1). 

The vector fields e 7 a are tangent to any r,t = constant two-sphere S, and so u a := e a 
and n a := e 1 a are, respectively, timelike and spacelike unit vectors orthogonal to S. From 
Eqs. ( [2.8|) we see that the orthonormal basis components of the extrinsic curvatures l ab and 
k ab are given by 

hj = ej a ej b V a u b = -u QJI , (5.8) 
ku = ej a ej b V a n b = -u UI . (5.9) 

Here ujqba = — e^°e B V a ec?6 are Ricci rotation coefficients. Working out these coefficients^ 
I find that the trace-free parts of Ijj and ku are given by 



where a = dgu, (5.10) 





ku = p(l °, ) , where (3 = ^d r In (5.11) 

Geometrically, the coefficient a is just e\ ■ f2(p recess ) ; i. e . ; the radial component of the angular 
velocity vector £}(p recess ) that measures the precession of a gyroscope carried by a locally non- 
rotating observer, relative to the observer's orthonormal frame (see Eq. (33.24) of Ref. p9 



compare also with Eq.( [3~2"l) above, and the discussion following it). Thus, the unreferenced 
shear term in Eq. ( |5.5|) is given by 

(P - P) = 2(/3 2 - a 2 ) = ^sin'e + O (1) , (5.12) 

where the last expression on the right hand side is the large r asymptotic expansion. So 
clearly, being of order 1/r 6 , (k 2 — I 2 ) does not contribute to the large sphere limit of the IQE 
at spatial infinity. 

What about the reference term {k — / 2 ) rcf ? It is plausible that the reference term is 
of the same order in 1/r as the unreferenced term, viz. 1/r 6 , or less, and so also does not 
contribute. However, to be certain one needs to solve the embedding equations ( |4.6H4.8j ), 
subject to the condition jF ref = JF, as argued in Sec. 4. We will not attempt to do so here, 
but it is instructive to at least work out what T is for the Kerr geometry. From Eq. ( |2.11| ) 



we see that the orthonormal basis components of the connection in the normal bundle are 
given by 

A! = ej b n c V b u c = -u m . (5.13) 

11 The easiest way to do this is to recognize that, with Z = or 1, ojzji — ctuj)Z: where a AB — i £B i eA de c . 
Thus we need only compute the exterior derivative of e 1 . 
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Evaluating these Ricci rotation coefficients reveals that the one-form A = Aje 1 (pulled back 
to the two-sphere S) is given by 



A = 1 d(f) } where 7 = ^Y r d r u = — sin 2 # + O — . (5.14) 

2Np r 2 \r 4 / 

Recall that u = u>(r, 0) is a measure of the frame dragging produced by the rotating geometry. 



While the shear in the time direction measures the 9 dependence of ui (see a in Eq. (|5.10 )) 



Eq. (|5.14 ) shows that the connection in the normal bundle measures its r dependence. Both 



are measures of angular momentum. 

Calculating the exterior derivative of A leads to the curvature in the normal bundle: 

T = - -d 6 — d r u = — cos 6> + — , (5.15) 



which is of order 1/r 4 . Inspection of the Ricci equation (|4.8| ), or Eq. ( |2.16|) with the left 



hand side set to zero, reveals that a solution to the reference embedding equations, subject 
to Eq. (|4.10|) , requires that l}j and kj e j be two matrices whose commutator is of order 



1/r 4 . On the other hand, one might guess that the trace of the difference of their squares, 
(k 2 — / 2 ) ref , might be of order 1/r 6 , as suggested above. It is not difficult to convince oneself 
that these two conditions are not incompatible, so the reference embedding equations at 
least do not obviously forbid the reference shear term (k 2 — / 2 ) rcf from being of the same 
order of magnitude as (k 2 — I 2 ), such that neither contributes to the IQE. 

In any case, assuming just that (k 2 — / 2 ) ref is at most 0<(r~ 3 ), which is almost certainly 
true, we find that the large sphere limit of the IQE at spatial infinity is given by 

lim IQE= lim — / dSraaR. (5.16) 

This limit of the IQE thus has a simple geometrical interpretation: apart from a factor 
proportional to the areal radius r, it is just the average over S of the sectional curvature of 
(S,a) as embedded in the physical spacetime (M,g). 

Now let us assume that (M, g) is vacuum (R a b = 0) near spatial infinity, so that there 
the Riemann tensor reduces to the Weyl tensor, C a b c d- From the definition of the two-surface 
metric given in Eq. fl2.3| ), and the fact that the Weyl tensor is traceless, one immediately 
gets 

aaR = 2E ab n a n b . (5.17) 

Thus the sectional curvature of (S, a) is just (twice) the radial-radial (Coulomb) component 
of the electric part of the Weyl tensor, E a b := —C ac bd,u c u d . Inserting this result into Eq. ( |5.16| ) 



we see that in this limit the IQE is precisely the coordinate-independent expression of the 
ADM mass given by Ashtekar and Hansen pU[.p| 

One more remark is in order here: Hayward's work on quasilocal energy || resembles 
what is done here in the sense that the Gauss embedding equation plays a central role, and 



12t 



In Ref. [39] a different definition of E a b is used, namely E a b = C ac bdn c n d ', but accounting for this difference 



in notation the two results agree. 
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that the analysis is boost invariant in spirit, i.e., no reference is made to a spacelike three- 
surface spanning S, with its attendant preferred time direction on S, and so on. However, 
Hayward's quasilocal energy is distinct from the IQE here: it does not involve a square 
root. Basically, Hay ward starts with an integral over S of the '2+2' Hamiltonian density, 
which yields a dimensionless quantity, and then multiplies this quantity by the areal radius 
of 5* to correct this 'defect', i.e., give it the dimensions of energy. This is in the same spirit 



as the areal radius factor appearing in the Hawking mass [ 25] . For the large sphere limit 
at spatial infinity Hay ward arrives at the same result given in Eq. (|5.16| ), except with r 
'outside the integral', so to speak. His quasilocal energy has the very appealing feature 
of not requiring a reference subtraction term, at least when the sectional curvature aaR 
falls off as 1/r 3 (however it diverges if, e.g., the spacetime is asymptotically anti-de Sitter 
space — recall Eq. ( |4.5|) ) ■ In our case, the square root in Eq. (|4.2j ) ensures that the IQE has 
the dimensions of energy, but the price paid is that a reference subtraction term is needed. 
(Without the square root the large sphere limit of the unreferenced IQE would just be 
(negative) the Euler number of S, which is finite, but carries no information about energy.) 
The areal radius factor r in Eq. (|5.16|) appears 'inside the integral': as mentioned before, it 



arises from the dominant scalar curvature term 1Z, and is really y2/1Z. I emphasize that 
this is a geometrically natural mechanism — r is not put in by hand. Finally, while one might 
feel that there is something unattractively ad hoc about a reference subtraction term, the 
flexibility it affords makes it possible to deal with the wide range of boundary conditions 
possible in general relativity. In Sec. 7 we will consider an interesting example of this. 

5.2 The null infinity limit 

We now suppose that the physical spacetime (M , g) is asymptotically flat at future null 
infinity. As in the previous subsection we begin with the generic large sphere form of the 
IQE given in Eq. (|5.5p, except now we take the large S limit in the future null direction. 



More precisely, on M we introduce the Bondi coordinates j40] x a = (w,r,8,(f)), and as 
before, denote the subset (8, 0) of spherical coordinates by x % . The retarded time w labels 
a one-parameter family of outgoing null hypersurfaces, and r is an areal radius (luminosity 
parameter) along the outgoing null geodesic generators of these hypersurfaces. The w,r = 
constant surfaces are topologically two-spheres, any one of which we denote as 5*. This setup 
is the same as discussed at the beginning of Sec. 5, where w here is what we there called the 
'time' coordinate, r. We are interested in the one-parameter family of two-spheres S in the 
limit as r — > oo, with w arbitrary but fixed. 

In the Bondi coordinates our asymptotically flat metric takes the standard form pT] . [T2[| 

g ab dx a dx b = -UVdw 2 -2Udwdr + a {j (dx i + W l dw) (dx j + W j dw). (5.18) 
We assume the following expansions for the various terms in this metric: 13 

V = l-2mr" 1 + 0<(r- 1 ), (5.19) 
We are following closely the notation used in Ref. p8[ , as well as the spirit of the discussion in their footnote 



13 



2. The meaning of the notation Q < (r n ) was described following Eq. ( |5.l| ). The motivation for this level 



of generality is that Chrusciel et al 42 have recently shown that one can allow 'polyhomogeneous' terms of 



r 



in these expansions and still have a consistent framework for solving the Bondi-Sachs-type 
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(5.20) 

(5.21) 
(5.22) 



The function V contains the mass aspect, m(w , 6 , <ft) . Observe that the metric on S is of 
the same form given in Eq. ( |5.1| ) (with Z = —X because r is an areal radius here), except 
now X(w, 6, <p) and Y(w, 8, 0) have a significant physical interpretation: they are the real 
and imaginary parts of Sachs' complex asymptotic shear c = X + iY ||41|| . Thus the scalar 
curvature of (S, a) will be given by Eq. ( |5.2|) , and we can begin our discussion of the IQE at 
Eq. (5J5). Our first task is to compute the unreferenced shear term, (k 2 — I 2 ). 



Inspecting the metric in Eq. ( |5.18|) , we choose the following basis of one- forms: 



e~ = ]-Udw, e + = dr + ]-Vdw, e 1 = j 1 i (dx i + W i dw) , 

_ 



(5.23) 



where indices I, J, . . . take the values 2 and 3, and , y I i is defined by demanding that cr^ 
\t 1 j- A suitable choice for 7^ is given by 





(5.24) 



In this matrix expression, I (i) is a row (column) index. Let the indices A,B,... take values 
in the set { — , +, J}. Then the metric is given by g a b = r]ABe A a e B b , where 



Vab 



( -200^ 

-2 

10 

V 1 J 



(5.25) 



This matrix, and its inverse, r\ AB , are used to raise and lower the the basis indices. The 
vector fields dual to the one-forms in Eq. ( |5.23| ) are given by e A a = rjABg ab e B b , or explicitly: 



2 
U 



d w - -Vd r 



W l di 



e+ = d r , ei = ^idi 



(5.26) 



dija^jJj, a l] being the inverse of the matrix Oi 



where 7/ is defined by 7/ — u 2J u - , •, w U uu 6 ux± c mvuoc ^ unc maunA u iy 

Inspection of e/ in Eq. ( ^.26| ) shows that these vectors are tangent to S, and so e± are two 
null normals to S. Since their normalization is such that e + ■ e_ = —2, we can set £± = e ± a , 

characteristic initial value problem. Allowing only expansions in powers of inverse r is tantamount to Sachs' 
'outgoing radiation condition' [ftl| , which they argue is overly restrictive. However, besides making the 
calculations 'tighter' as regards remainder terms, and slightly more general, we would get the same results 
had we assumed Sachs' outgoing radiation condition. 
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where £± was previously defined by £± := u a ± n a (see Eq. fl3.8Q ). Thus, from Eqs. ( |2.8j ) we 
have the following result, in basis components: 

lu ± kjj = ej a ej b V a e± b = -u ±JI , (5.27) 

Working out the required Ricci rotation coefficients^ I find 

llj + k U = B U , llJ -k U =Tr ( A U ~ " V d W J)) > ( 5 - 28 ) 



where 



U 



Au = 7 ( /7^ = ; ( f ^ ) + ° <(r_1) ' (5 ' 29) 

Here we use an overdot (prime) to denote differentiation with respect to w (r). Observe that 
T>nWj\ in Eq. ( |5.28 ) is of order 1/r 2 . Taking the trace-free part of Eqs. Q5.28D gives us the 
basis components of the shears in the two null directions: s±u = lu ± kij- Explicitly: 



s+u 



X Y \ „ , 2 (X Y 



r 



0<(r~ 2 ), s_u = - I V )+0<(r- 1 ). (5.31) 



2 v y -x / 1 ^ v 7 ' JJ r \ y -x 



Thus the unreferenced shear term in Eq. ( |5.5|) is given by 

(P - f 2 ) = -s +IJ s_ IJ = ^(XX + YY) + 0<(r- 3 ) (5.32) 

(cf. Eq O). 

We thus learn that, in contrast to the spatial infinity limit (see Eq. ( j5.12j )), in the null 
infinity limit the unreferenced shear term is of order 1/r 3 , and so does contribute to the IQE. 
We will argue below that the reference shear term (k 2 — l 2 ) rel is also of order 1/r 3 , but that 
it is a total derivative and therefore does not contribute. So as not to interrupt the flow or 
our discussion, for the moment let us assume this is true, in which case Eq. (|5.5|) becomes 



lim IQE = lim / dS r 

r^oo r^oo lQ n Jo 



ooR - 4r(XX + YY) 



(5.33) 



Because of the cc + cc = 2 (XX + YY) term, this result looks like it could very well be the 
Bondi-Sachs mass [[|l|] . To see that in fact it is, a straightforward calculation of the Riemann 
tensor of g a b projected into S gives the following sectional curvature of (S, a): 

ooR= -|(m + XX + YY) + 0<(r- 3 ). (5.34) 



See the footnote on page E3l with Z = ±. 
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Inserting this result into Eq. ( |5.33j ) we see that the shear terms cancel, leaving only the mass 
aspect, m: 

limIQE = — / dflrn(w,e,(b). (5.35) 

r->oo 4tt Js 

In obtaining this expression, recall that because r is an areal radius we can (and did) take 
yfa = r 2 sin 9, and so dS = r 2 dQ, where dfl = sin 9 d9 dip is the measure on the unit round 
sphere.^ Thus the future null infinity limit of the IQE is the Bondi-Sachs mass |41 . 



Now there is an important lesson to be learned from this result. The unreferenced shear 
term (k 2 — I 2 ) is solely responsible for producing the all-important cc + cc term that accounts 
for the mass loss due to gravitational radiation. Hence this term is necessary under the 
square root in Eq. ( [4.2D , and so there is no natural way to avoid (k 2 — / 2 ) ref in IQE ref , and 
its attendant embedding problem. Moreover, we learn that these shear terms are not only 
associated with angular momentum, as I have been stressing, but also encode information 
about gravitational radiation. We will see precisely the same phenomenon emerge in the 
small sphere limit in Sec. 6. Furthermore, it is emphasized in Ref. [[|3|] that it is easy to 
construct, ab initio, an integral expression involving the Riemann tensor (e.g., an integral of 
ao~R over S) that is conserved under certain circumstances. One is thus tempted to interpret 
such a conserved quantity as an energy. However, such attempts fail to produce, in the null 
infinity limit, the crucial null- surface- dependent shear terms seen in Eq. (|5.33| ), and it is 
difficult to see how to modify them in a covariant way to produce these terms [J43|J . The 
shear term (k 2 — I 2 ) is precisely such a covariant modification. Moreover, it arises naturally 
from simply replacing the Brown- York k with the boost invariant quantity \Jk 2 — I 2 . (Of 



course a similar observation can be made concerning, say, the Hawking mass |25| , which has 
the same large sphere limit as in Eq. ( [5.5|) , except without the reference shear term.) 

These clean results rely on our assumption that the reference shear term (k 2 — / 2 ) ref does 
not contribute to the null infinity limit of the IQE. I claimed above that this is so because it 
is a total derivative. To prove this would require solving the embedding equations ( [4.6H4.8j ), 
which we know is a very difficult task. However, I will now present a heuristic solution of 
the Ricci equation that leads to a substantiation of this claim. Moreover, we will see how 
demanding = T plays a crucial role in achieving this result, which provides our first 
bit of indirect but concrete evidence that this condition is required to properly account for 
angular momentum (and as we now see, also gravitational radiation). 

To begin we need to calculate the connection in the normal bundle, A a , and then its 
corresponding curvature, T '. Proceeding as we did in the spatial infinity case (see Eq. ( |5.13| )) 
we find that the basis components of A a are given by 

Aj = ±u + -j = ^jWj + i ej (lnZ7), (5.36) 

where ej(lnU) denotes the derivative of InU along the vector field ej. This term is pure 
gauge. As for the other term, since we only know Wi to leading order we can put U = 1 



15 



On this note, it might be helpful to point out an important detail in the calculation of the sectional 



curv ature given in Eq. ( 5.34 ). One of the the terms that arises in the calculation is the trace of Ajj in 
Eq. ( pT29|) , which appears to be of order < (r 1 ). If this were so it would be problematic. But in fact it is 



zero, because A/ = r y I ' lA f I i = (V^)^/ " = 0, since a = r 4 sin 2 6 does not depend on w. 
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here — see Eqs. ( p.2U| ) and ( [5.211 ). Thus, up to a gauge transformation, the connection one- 



form A is just (1/r times) the one-form Widx 1 . And the curvature is thus proportional to 
the curl of W: 

T = ^e IJ VjWj = -(V 2 W 3 - V 3 W 2 ). (5.37) 

Keep in mind that the numerical indices here refer to basis components, not coordinate 
components. It is easy to see that T is of order 1/r 3 . It is interesting to compare T with 
the scalar curvature 1Z, whose form was given in Eq. (|]D. Using the metric in Eq. (|5T22p it 
is not difficult to evaluate the remainder term, An- The net result is 



2 2 



K = ^ + -P-iy + 0<(r- 3 ), (5.38) 



r 



Note from Eq. (|5.21 ) that W i is of order 1/r 2 , so the term above involving W is, indeed, 



of order 1/r 3 , as it should be. Thus we see that TZ is associated with the divergence of W, 
and T with its curl. This is an explicit example of a point made earlier, namely that both 
curvatures are on the same geometrical footing: To capture the two pieces of information in 
W — its divergence and its curl — requires precisely both 7Z and T . 

Let us now turn to the null shears of S as they appear in the physical spacetime, 
Eq. ( [5.31| ). The form of these shears suggests we make the following ansatz for the null 
shears of S in the reference (Minkowski) spacetime: 




^ref I P I I {ry,~2\ n re f 




+ 0<(r" 2 ), = - I " +0 < (r- 1 ). (5.39) 



Observe that we might expect the pair (a, (3) to play a role distinct from the pair (7,(5). 
Comparing Eqs. ( |5.39|) and ( [5.31 ) suggests that a and f3 will be like X and Y in that they 



have something to do with the intrinsic geometry of (S, a). The 'more important' terms will 
be 7 and 5, because they occur at the dominant power of inverse r. Also, we expect them to 
be related to the extrinsic geometry of S, since their counterparts, X and Y, measure how 
aij changes as a function of the retarded time w — they are the two 'news' functions |[41|| . 

With these observations in mind, we will now 'solve' the Ricci embedding equation, (|4.8| ). 
We first compute (in basis components) 

" r'jK'i sr -Vi = - M) = ^ 3 det ( a P ) . (5.40) 



Next we impose the condition jF rcf = JF, and observe that T in Eq. (|5.37| ) can also be 
expressed 'determinant', i.e., 

>" = w 3 ) (5 ' 41) 

The Ricci embedding equation instructs us to equate the two previous determinant expres- 
sions. One solution is to make the identifications a «-> rT>2 and (3 <-> rT> 3 (which are 
consistent with our expectation that a and (3 be associated with intrinsic geometry), to- 
gether with 7 <-> rW 2 and 5 <-> rW% (which are consistent with 7 and 5 being associated 
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with extrinsic geometry, since W is proportional to the connection in the normal bundle — a 
measure of extrinsic geometry). Notice that this means s+jj is a derivative operator! To 
make this more palatable one may go to a 'Fourier transform space', where the derivative 
operators T>j become momenta, /C/. Accepting these heuristic identifications, and recalling 
Eq. (|4.9|) , the key point now is to observe that 

(P _ J2yei = ref / ref J = + fig) = -1(V 2 W 2 + V 3 W 3 ) = --V • W, (5.42) 

the result we desired: being a total derivative, [k 2 — Z 2 ) ref does not contribute to the IQE. 

Now recall that the purpose of solving the embedding equations is to establish a rela- 
tionship between jF rcf and (k — Z 2 ) ref . It is via this relationship that we envision the angular 
momentum information in jF rcf (= JF) to enter IQE ref , which sees only (k 2 — Z 2 ) rcf . The 
previous result implies the following relationship: is the curl of W, and (k 2 — Z 2 ) ref is the 
divergence of W. Hence, no relationship! Is this a problem? Certainly not. We require the 
relationship in question only when different isometric embeddings of the same (5*, a) result 
in a different {k 2 — Z 2 ) ref , making the IQE ambiguous. Such a relationship can be used to 
resolve this ambiguity by setting jF ref = T . But accepting the argument that ik 2 — / 2 ) ref 
is always a total derivative, any such ambiguity would be of no consequence here. Further- 
more, in this simple case there is no need for ik 2 — Z 2 ) rof to carry any information about 
angular momentum (or gravitational radiation) coming from J 7 , because all of the relevant 
information is already carried in the unreferenced shear term, (k — I 2 ). This is not to say 
that the condition jF rcf = T is not important here. According to our heuristic argument, it 
is precisely this condition which ensures that (k 2 — Z 2 ) ref is a total derivative, and so carries 
no information. 

The null infinity limit is a simple case which only minimally exercises the consequences 
of the condition jF ref = T. For (S, a) a finite two-surface the situation is much richer: it 
is highly nonlinear, since the square roots do not disappear, and almost certainly requires 
effectively a one-to-one relationship between T and (k 2 — l 2 ) rei to remove the embedding 
ambiguity inherent in the IQE. 

6 The small sphere limit of the IQE 

Having considered the large sphere case, we now turn our attention to evaluating the 
IQE when (S, a) is a small sphere. The large and small sphere limits are similar in that 
in both cases S approaches an asymptotically flat region of (M,g). In the latter case, the 
asymptotically flat region is the infinitesimal neighborhood of a generic spacetime point 
p G M, which is the 'center' of our shrinking sphere. For simplicity we will suppose that 
(S, cr) is asymptotically round. Another feature in common with the large sphere limit is 
that in this codimension-two setting, the limit can be approached from different directions, 
either spatial or null. More precisely, fix a set of Riemann normal coordinates (t, x' 1 ) about 
the point p, set r 2 := 5ijX l x^, and define S 1 * by the condition (r, t) = (r*, ar*), where a is a 
direction parameter. Then consider the limit of S 1 * as r* — > 0. As before, we will henceforth 
omit the subscript '*'. The case a = is a spatial limit, since then S always lies entirely in 
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the t = spacelike three-surface containing p. a = ±1 is the null limit, in which S lies in 
the future/past light cone of the point p. The latter case was considered by Horowitz and 
Schmidt p0| in their classic work on the small sphere limit of the Hawking mass. Brown, 



Lau, and York [TJj also consider this same limit of the Brown- York quasilocal energy. We 
will be borrowing some results from these two references. 

Explicitly, for a given value of the parameter r, S is defined as a submanifold of (M,g) 
by embedding a topological two-sphere with coordinates 9 and into the Riemann normal 
coordinate system, as follows: t = ar, x 1 = r sin 9 cos <f), x 2 = rsin^sinc/), x 3 = rcos9. Since 
(with t = ar) the deviation from the flat metric in Riemann normal coordinates is 0(r 2 ), 
the induced metric a a b on S will differ from that of the round sphere to this same order, and 
so the scalar curvature of (S,a) will have an expansion in r of the form 

K = \ + ft<°> + rll^ + r 2 n^ + 0(r 3 ), (6.1) 

where each of the coefficients is a function of 9, 0, and the parameter a. To evaluate 
the IQE we will also need similar expansions for the other quantities appearing in Eqs. ( |4.2[ ) 
and ( |4.3|) . These are written as follows: 

aaR = aaR {0) + r aaR w + r 2 aaR {2) + 0(r 3 ), (6.2) 
(P-[ 2 ) = r 2 (P-/~ 2 )^ + 0(r 3 ), (6.3) 
(P _ p) = r 2 {k 2 - l 2 f^ ci + 0(r 3 ), (6.4) 

where each of the coefficients on the right hand side is similarly a function of 9, <f), and 
a. Since the appropriate reference spacetime in this case is Minkowski space, we have 
aaR rei = 0. Substituting these expansions into Eq. (|4.2j) we find that in the small sphere 
limit the unreferenced IQE behaves as 

IQE unref = _ J_ J '^dS^ll + j [(TZ W - aaR {0) ) + r(K {1) - aaR m ) (6.5) 



+r 2 ( n (H) _ aaR (2) + ( p _ p )( 2) _ I ( ^(0) _ aajR (0) )2 



+ 0(r 



5 a 



Similarly, the reference IQE behaves as 
1 f ,„2|\ r 2 

?7T JS 



8tt Js r{ 4 L 

+r 2 (n^ + (P - f»)(2)«* - ^(7^) 2 )] + 0(r 5 )| . (6.6) 

Notice that, unlike in the large sphere limit, neither the unreferenced nor reference energies 
diverge as r — > 0. Nevertheless, the reference subtraction procedure is still necessary to 
eliminate the leading 0(r) term in IQE unref , which has nothing to do with energy. Thus, 



31 



forming the difference of the previous two expressions we find that the small sphere behavior 
of the (referenced) IQE is given by 



IQE 



167T JS 



dSr 



aoR - (P - I 2 ) + (P - P) rcf + -r 2 aaR^\aaR {0) - 27Z {0) ) + 0(r 3 



(6.7) 
.2W. 



The sectional curvature and shear terms have been resummed according to Eqs. 
and the resulting expression is valid to the order indicated. Notice that, as in the large sphere 
limit, the scalar curvature TZ dominates the other terms under the square root, allowing us to 



expand the radical about y4/r 2 = 2/r.0 And after the reference subtraction is performed 
what remains again is the sectional curvature of (S, a) as the dominant term contributing to 
the energy. Comparing Eqs. ( |6.7]) and ( |5.5| ) we observe that both the small and large sphere 
limits of the IQE are very nearly formally identical. 

We will split the IQE into three pieces, each to be discussed separately: IQE = IQE X + 
IQE 2 + IQE 3 + 0(r 6 ), where 



(6.8) 
(6.9) 
(6.10) 



IQE 1 = 


16tt j 


f dSr 

s 


aaR - (P 


- I 2 )] , 


IQE 2 = 


1 

1287T 


[ dSr 
Js 


3 [aaR^ \aaR^ -27Z^)} 


IQE 3 = 


16tt j 


[ dSr 

s 


"(P - i 2 y cf 





We begin with IQE 1; and show that this piece is essentially the Hawking mass P5| . To see 
this, we combine Eqs. ( fj~T|) and ( [2.10|) to get 



ooR - (P - I 2 ) =TZ- -(P - I 2 ) = TZ - 2H ■ H. 



(6.11) 



Replacing the integrand of IQE 1 with the last expression, and using the Gauss-Bonnet the- 
orem to integrate the TZ term, we find 



I QEl= * M 
^ 1 4tt V 4vr 



2ti- - f dSH-H 
2 Js 



(6.12) 



where we pulled r outside the integral and replaced it with A/(4tt), where A is the area of 
(S, a). This form of IQE : is precisely the expression of the Hawking mass given in Ref. pOf . 
(Our H° in Eq. ( [2.10|) is their N c /2, and the sign of their metric signature is opposite to ours.) 
Comparing Eq. ( |6.12[ ) with Eqs. ( |3.7]) and ( |3.9| ) we observe that, while the unreferenced IQE 
involves the mean curvature itself, y/H • H, the Hawking mass is constructed from the square 
of the mean curvature. As mentioned above, the square root in y/H ■ H effectively disappears 
in the small (and large) sphere limits due to the presence of the dominant scalar curvature 

16 As remarked in the footnote on page if the sphere is not asymptotically round the shear terms will 
contribute at order 1/r 2 , and the 2/r term outside the braces in Eqs. ( |6.5| ) and ( |6.6| ) will be modified 
accordingly. We will not consider this more complicated case here. 
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term, and consequently the leading order contribution to the IQE reduces to essentially the 
Hawking mass. 

There are two subtleties worth mentioning: (i) Replacing r with ^ A/(An) is in general 
not valid because it requires that r be an areal radius which, in general, it is not. However, 
it certainly is to lowest order in r, which will be sufficient for our purposes here. But to 
higher order, IQE X and the Hawking mass will in general give different results, (ii) It is well 



known that the Hawking mass runs into difficulties when (S, er) is not a round sphere [H,E0 



a problem that was addressed by Hayward in Ref. ||. It might be that this problem is a 
result of having to insert by hand the factor \J A/(Aix) outside the integral, versus having r 



'inside the integral' generated automatically by J 2/11. A related issue was discussed at the 
end of Sec. 5.1 in connection with Hayward's definition of quasilocal energy. 

The connection between IQE : and the Hawking mass allows us to borrow some results 



from Ref. [2"0|, which are calculated for the null limit case (a = 1). When matter is present 



Horowitz and Schmidt find that, to lowest order in r, the Hawking mass is 

IQ E i = (Jvrr 3 ) TTu a u% + 0(r 4 ). (6.13) 

Here T^ at is the energy-momentum tensor of matter, and the expression is to be evaluated 
at the point p, where the unit timelike vector u a is just (d/dt) a in our Riemann normal 
coordinates. This is a standard result in the literature on quasilocal energy [pj|,|20|,|2l|. [22|. 



23, 2"3}|, and a very significant one. As emphasized in the Introduction, the quasilocal idea 
asserts that the time-time component of the energy-momentum tensor of matter a priori has 
nothing to do with energy. It is only from the small sphere limit of the quasilocal energy that 
we learn this quantity is an energy volume density, i.e., multiplying it by the volume factor 
47rr 3 /3 gives the energy in an infinitesimal sphere of proper radius r. However, integrating 
this energy volume density over a finite volume to determine the total energy inside is not, 
in general, valid unless one wishes to ignore gravitational effects which, as we will see in 
moment, come at higher order in rf£] It is in this sense that the quasilocal idea implies that 
even energy due to matter is not localizable in the context of general relativity. 

Now let us assume the spacetime is vacuum in the neighborhood of p. Then the leading 
order contribution to the Hawking mass is |2(| 



IQEj = ^T abcd u a u b u c u% + 0(r 6 ), (6.14) 



where T abcd is the Bel- Robinson tensor g5j. Thus gravitational energy begins to appear at 
0(r 5 ). This same result is obtained for the Brown- York quasilocal energy for a suitable choice 



of reference embedding [[L9j. However, this is not a universal result in the literature [21.22 



23|,|24[] . For example, Hayward's quasilocal mass gives a similar result as above, but with the 



17 A well established example of this phenomenon is the Tolman density, which integrates to the Komar mass, 
and is defined in the special case that the spacetime is stationary and asymptotically flat jjj . It is noteworthy 
that it is not T a b that appears in the Tolman density, but rather the combination T a t — (l/2)Tg a i,. The extra 
term involving the trace of T a f, is associated with gravitational effects — see Ref. Q| , and problems 4 and 5 
in Chapter 11 of Ref. §J. 



33 



numerical factor 1/90 replaced with —2/45 Given that gravitational energy is such a 
difficult problem it is not surprising that a consensus has not yet been reached. 

We now turn our attention to the second contribution to the IQE, namely IQE 2 given 
in Eq. ( |6.9| ). This quantity represents a deviation from the Hawking mass due to the fact 
that the IQE is roughly the square root of the former. Actually, it is clearer to compare 
with the Brown- York quasilocal energy, since in constructing the IQE we simply replaced 
the Brown- York k with \Jk 2 — P. In the context of our generalization given in Eq. (|3.4j ), we 
therefore have the heuristic comparison: 

m = JE 2 -p 2 = E- — ► — Vk 2 -P = — (k- — V (6.15) 

v y 2E 8tt 8tt V 2k J K ' 

So IQE 2 might be thought of as the analogue of the term —I 2 /(2k), and as such would be 
expected to reduce the magnitude of the IQE from the result given in Eq. ( ft.14 ). 



In order to calculate IQE 2 we need to evaluate the quantities and aaR^ in Eq. (|6T 
To do so we appeal to the Gauss equation, (|4.1|). Up to zeroth order in r, this equation reads 



oaR^ = - + - -(k 2 - I 2 ) + Otr), (6.16) 
r 2 2 



where we made use of Eq. (|6.1|) . We will show later that 



4 4 
3' 



(k 2 - I 2 ) = — + -(1 + 2a 2 )E ab n a n h + O(r), (6.17) 



where a is the direction parameter introduced at the beginning of this section. E ab n a n b is 
the radial-radial component of the electric part of the Weyl tensor, which we saw earlier 
in Eq. ( |5.17|) . This quantity (E a bn a n b ) is to be evaluated at the point p, where in our 



Riemann normal coordinates the radial unit vector n a has only spatial components, given 
by n l = x l jr. We also need dS to lowest order, which is just r 2 dQ, dfl being the measure 
on the unit sphere. Putting these results together we have 

IQE 2 = —(5 + Aa^EijEki f dfinWnV. (6.18) 

967T ' J 

By symmetry, the integral over the product of radial vectors must be proportional to 5^S kl \ 
Transvecting both this term and the integral in question with S^Sm, we easily obtain that 
the proportionality constant is 47r/5. Using the fact that E a b is symmetric, trace- free, and 
orthogonal to u a , we find that 

IQE 2 = --i-(5 + 4a 2 yE ab E ab . (6.19) 
loll 

So IQE 2 is negative, as expected, and this negative contribution is to be added to IQE : in 
Eq. (|6.14|) . Of course we can only consider the a = 1 case, since this is the case assumed in 
Eq. (|6.14j) . Recall that the time component of the Bel- Robinson tensor can be expressed in 



terms of the electric and magnetic parts of the Weyl tensor []45 |: T a b c dU a u u c u = E ab E + 
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B a bB ab , so IQE X is nonnegative. Inspection of Eq. ( |6.19| ) shows that adding to IQE X the a = 1 
value of IQE 2 makes the energy have indefinite sign. It is positive (negative) if the magnetic 
(electric) part dominates. This seems like a strange result, but it is only an intermediate 
result. We have not yet considered the last contribution, IQE 3 , involving the reference shear 
term. But unfortunately at present I do not know how to solve the embedding equations to 
determine this term. 

Now one can construct a heuristic argument much like the one given at the end of Sec. 5.2, 
which suggests that {k 2 — Z 2 ) ref is a total derivative, and so does not contribute. However, 
the argument is much less believable in this case. In contrast to Eq. ( |5.39| ) it turns out that, 
because T is 0(1) in r (as we shall see later), we must expand the reference shears s± over 
three orders of magnitude in r, from 0(1) to 0(r 2 ). One might trust a heuristic argument 
working to leading order, but believing the higher order corrections is less palatable. In 
short, I do not know what energy prediction the IQE gives at 0(r 5 ), and until a solution to 
the embedding equations is found there is no sense in speculating. 

However, before leaving this section I will provide an intriguing interpretation of how a 
definition of quasilocal energy such as the Hawking mass (or the IQE) provides a measure 
of the gravitational energy contained inside a small sphere. 

In the a = 1 null limit case, the lowest order contribution to the gravitational energy, 
namely (l/90)r 5 T a b c< iU a u b u c u d \p in Eq. ( |6.14| ), originates in the 0(r 2 ) terms inside the brack- 



ets of the integrand of IQE X in Eq. ( |6.8| ). In the terminology of Eqs. ( |6.2|) and ( |6.3| ), this 
means we are interested in the coefficients aaR^ and (k 2 — l 2 p 2 '. Inspecting the Appendix 
of Ref. [19] reveals that these two coefficients differ only by a numerical factor. They are 
both proportional to ^o^olp (i n Newman- Penrose notation), and the two numerical factors 
conspire to produce the 1/90 factor in the final result. Thus, to understand how the inte- 
grand of IQE X encodes information about gravitational energy it suffices to study the shear 
term, {k 2 - l 2 )^ . We will now compute this term for arbitrary a G [—1,1] to see how it 
behaves in both the spatial and null limit cases of the small sphere. 

Denoting our Riemann normal coordinates {t,x l ) collectively as x a , the metric in these 
coordinates take the form 



9ab = Vab ~ \ > JabcdX C X d + 0(x 3 ), 



(6.20) 



where J a bai = (Racbd + Radbc)/2 is the Jacobi curvature tensor P9]| . We first construct a 
pair of mutually orthogonal unit normal vector fields u a and n a on S, with u a normal to the 
t = constant surface passing through S. These are given by 



u =-, u 



--r 2 /?° + 0(r 3 ), Uo = -N, = 0, 



(6.21) 



n 



0, n l = p 



,-vt ^ j ry 

^ + V/3^ + 0(r 3 ) 
r 3 r 



no = ~r 2 p 0j - + 0(r 3 ), n % = p-, (6.22) 

O T T 



where 



N = l + -r 2 (3 00 + O(r 3 ) 
6 



(6.23) 
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P=l-V/%^ + 0(r 3 ), (6.24) 
6 r z 

Pab = a JabOO + 2aJ a fc j— + J a bij—^~- (6.25) 

Since the Jacobi tensor in Eq (|6.20| ) is evaluated at the coordinate origin p, its indices, 
and thus those of P a b, are raised and lowered with the flat spacetime metric r]^ = rf h = 
diagonal(— 1, 1, 1, 1). Similarly, Xi := S^xK 

Now define on S a pair of mutually orthogonal unit tangent vector fields e 7 a , where 
indices I, J, . . . take the values 2 and 3. The set {e a := u a , e x a := n a , ej a } thus comprises 
an orthonormal basis adapted to 5*. Let basis indices A,B,... run from to 3, and a, /3, . . . 
from 1 to 3. Beginning with this setup it is straightforward to compute the basis components 
of the extrinsic curvatures defined in Eqs. Q2.8|) . I find: 

l u = ej a e^ a u h = ~r[aJ 00I j + J mIJ ] + 0(r 2 ), (6.26) 

' ' ' ' Xl + 0(r 2 ). (6.27) 



ku = e/e/V n 6 = ~5 U - -r 
r 3 



In these equations a quantity such as Johj means [e a e l b e I c e J d Jabcd\\p, which is a function of 
only the angles 9 and <fi on S. 

Since we are interested in purely gravitational energy we shall restrict ourselves to the 
vacuum case. In our basis components the electric and magnetic parts of the Weyl tensor 
are defined by 

E a /3 = —Co a oi3 and B a p = — *Co a o/3, (6.28) 

where *Cabcd = {^I^)^ab EF ^efcd- These are symmetric trace-free three-dimensional 
tensors associated with t = constant spacelike hypersurfaces. As such, each has five inde- 
pendent components, which together comprise the ten independent components of the Weyl 
tensor. In terms of these fields, the components of the Jacobi curvature tensor relevant to 
Eqs. ( |6T26| - |6T27D read: 

-Toon = — En, Joou = —Eu, Joiu = —*B IJl Jmj = —Ejj — EuSij. (6.29) 

Here Bu is the trace- free part of Bjj, and *Bu = ej K Bxj is its dual in (S, a), which is 
also trace- free. The trace of Eu is 5 IJ Eu = b a ^E a p — En = —En, since E a p is trace-free. 
Thus, the trace of the extrinsic curvatures in Eqs. ( |6.26| - |6.27| ) is found to be 

I = -~raE n + 0(r 2 ), (6.30) 
3 

k = - + \r{l + 2a 2 )En + 0{r 2 ). (6.31) 
r 3 

Squaring these and forming their difference leads to Eq. (|6.17| ) written earlier. The trace-free 
parts are 

lu = \r{aE u + *B U ) + 0(r 2 ), (6.32) 

~kij = \r{l-a 2 )E L j + 0{r 2 ). (6.33) 
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Now Eu and Bjj each have two independent components, and it is useful to introduce the 
notation 



\{B 22 — E 33 ) E 23 \ f e 



.2 e s 



Eu= 2V " MJ \ J=: Is „2 , (6-34) 



£"23 —2 (-E22 — £33) / \ e 3 — e 

s _( \{B 22 -B 33 ) B 23 \ /6 2 6 3 \ 

5/J " I B23 ~\{B 22 - B 33 ) ) - ~ [b* -b 2 )> ^ 

and thus define a pair of two-vectors e := e I e 1 a and b := 6 / e / a tangent to S. Since e and 
6 represent the pullbacks to S of E a/3 and £? a/ g, respectively, they are to be thought of as 
'electric and magnetic fields' induced on S by the Weyl curvature that S is embedded in. It is 
easy to see that under a rotation of the basis vectors e 7 a through an angle 7, the components 
of e and b rotate through an angle 27, so e and b are not true (spin one) vectors, but rather 
spin two objects, as one would expect. 

Thus we arrive at the results we are interested in. To lowest order in r we have 

P = \r 2 {aE + *B) 2 = \r 2 (a 2 e- e + b-b-2aex b), (6.36) 
9 9 

~k 2 = \r 2 {l - a 2 ) 2 E 2 = \r 2 {\ - a 2 fe- e. (6.37) 
9 9 

The difference of these two gives the 0(r 2 ) piece of the unreferenced shear term appearing 
in the integrand of IQE X in Eq. fl6.8| ). Notice that it is the appearance of *Bjj (rather than 
Bij) that gives rise to the cross product term ex 6 = e 2 b 3 — e 3 b 2 in Eq. ( |6.36| ). 

We first consider the case a = 1, in which S lies in the future light cone of the point p. 
Then k 2 = and so the shear term (k 2 — I 2 ) is proportional to r 2 (E + *B) 2 . As a quick check, 
it is easy to verify that (E + *B) 2 , in turn, is proportional to ^o^olp? i n agreement with 
Ref. |19[]. I mentioned above that in this case the aaR term in Eq. (|6.8|) also contributes a 
term proportional to r 2 \E f ^ r o|p |^|- Putting in the numerical factors I find that 



r r 3 r 1 / -> 1 



,e ■ e + b-b) e x b 

8vr v ; 4tt 



+ 0(r 6 ). (6.38) 



Now (e ■ e + b ■ b)/(8ir) looks like the energy surface density of the 'electromagnetic field', 
but we must be careful about its dimension, (e • e + b ■ 6)/(87r) has dimension L~ 4 , where L 
means 'length', which is not correct. However, the additional factor of r 3 /9 in the integrand 
suggests that it is really £ := r 3 (e • e + b ■ 6)/(727r) that is the proper energy density. £ 
has dimension L" 1 , consistent with it being interpreted as the 'electromagnetic' energy per 
unit area of S. Besides giving the right dimension, the additional r 3 factor forces f s dS £ 
to go to zero as r 5 , consistent with the fact that there can be no gravitational energy at 
order r 3 . We interpret J s dS £ as the total 'electromagnetic' energy that was in S at t = 0, 
or equivalently, the total gravitational energy that was in the small volume spanning S at 
t = 0. 

To further justify this interpretation we now turn to the radiation term in Eq. ( |6.38| ). 
Clearly ex b/(4ir) might be thought of as the gravitational analogue of the electromagnetic 
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Poynting flux, directed radially outward from S. But again, its dimension is wrong. Of 
course the factor of r 3 /9 will fix this problem, as before, but the situation is more interesting 
this time. Multiplying by r 2 /9 we get the proper Poynting flux, V := r 2 e x 6/(367r). V 
has dimension L~ 2 , consistent with interpreting it as the 'electromagnetic' energy per unit 
time per unit area. So J s dSV gives the 'electromagnetic' energy per unit time radiating 
from (or through) the surface S. The factor of r 2 indicates that the efficiency of a small 
volume to radiate gravitationally grows in proportion to its surface area, in analogy with 
an electromagnetic antenna. But there is one more factor of r, which one might imagine is 



the r outside the brackets in Eq. ( |6.8|) , i.e., y2/1Z. This distinction between rs is suggested 
by the close analogy between V and the shear term responsible for radiation in the null 
infinity limit — see Eqs. (|5.32[ - |5733|) . This additional factor of r has the interpretation of a 
time lapse, i.e., r J s dSV is the amount of 'electromagnetic' energy radiated from S between 
time t — and t = r. Thus, the following picture has emerged regarding Eq. ( |6.38| ). 
The 'electromagnetic' energy in S (or equivalently, the gravitational energy in the volume 
spanning S) at time t = r is the energy at t — minus the amount of energy radiated during 
this time interval. (Keep in mind that e and b are evaluated at p, and hence at t — 0.) 

The case a = — 1 is similar, except now S lies in the past light cone of the point p. 
Inspection of Eq. ( 16.36 ) reveals that the radiation term in Eq. (|6.38 ) now appears with the 
opposite sign. The fact that this sign change comes out correctly is reassurance that our 
picture is correct: The energy at time t = —r is the energy at t — plus the amount of 
energy that will be radiated from the sphere during the interval from t = — r to t — 0.Q 

Finally, we consider the spatial limit case, a = 0. According to our discussion above we 
would expect IQE : to be the same as in Eq. ( |6.38 ), except with the radiation term absent. 
Inspection of Eqs. ( |6.36[ - |637D reveals that this is not the case. However, it is only when a = 1 
(and presumably also when a = — 1) that we know that aaR^ in Eq. (|6.8j ) is proportional 
to (k 2 — l 2 )( 2 \ in which case it suffices to consider only the shear term. Unfortunately, it is 
not possible to compute aaR to 0(r 2 ) within the framework of our 0(r 2 ) Riemann normal 
coordinates, so we cannot learn if this simple proportionality between the two persists when 
| at | < 1. One might guess that it almost certainly does not, but I will leave this question for 
future work. Nevertheless, since we expect I 2 to play the key role with regards to radiation, 
Eq. ( |6.36|) is still of some qualitative value when |at| < 1. From this equation we see that the 
radiation term is zero when a is zero, and turns on in proportion to a, precisely as it should 
since the time lapse is now or, instead of r. 

As satisfying as it might seem, the picture just given is not truly quasilocal, in the sense 
that e and b are evaluated at the point p. To be truly quasilocal we need 'electromagnetic' 
fields, call them E and B, evaluated on S. This is where observers reside, and measure- 
ments are made, according to the quasilocal idea. Such a quasilocal picture is achieved very 

18 The reader may have noticed that b in Eq. ( |6.35| ) was defined with an awkward minus sign. This sign was 
chosen to give the picture just described. Reversing the sign is equivalent to replacing a with —a. Insofar 
as e and b (like E and B in electromagnetism) are defined by their physical interpretation, choosing the sign 
of b to give a result with the correct interpretation is legitimate. But this assumes we know what the correct 
interpretation is, and it is not certain we do. For example, I mentioned above that at 0(r 5 ) Hayward's 
quasilocal energy gives a negative gravitational energy |22|. If this is correct, then we should replace the 
definition of b with —b. 
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naturally as follows. The basic idea is: e and b are certain components of the Weyl tensor 
evaluated at p (r = 0). But this information is contained in the 0(r) piece of certain con- 
nection coefficients evaluated on S (r > 0). Thus we expect the desired E and B fields to 
be associated with connection coefficients. 

For simplicity we will restrict ourselves to the case a = 1, which also allows us to borrow 
some results from Ref. |Tj|. We first observe that 

-2[(e 2 -b 3 )+t(e 3 + b 2 )}. 



0\p 



(6.39) 



On the left hand side is a component of the Weyl tensor in Newman-Penrose (NP) notation. 
Using Eqs. ( |6.28| ) and ( |6.34j - |6.35| ), \&o is easily converted to the expression given on the right 
hand side. From Eq. (B5b) of Ref. |TJ| we have 



a = -^ | P + O(r 2 ), 
where a is one of the NP spin coefficients. Thus we find that 



(6.40) 



-aa = — [e • e + b ■ b - 2e x b] + 0(r 3 ). 
Comparing this with the integrand in Eq. (|6.38| ) we are led to define 



E := -e + 0(r 2 ) 
3 



B := -b 
3 



0(r 2 



or in other words, 



a = -2[(£ 2 - B 3 ) + i{E z + B 2 )] + 0{r 2 



Thus, to 0(r), E and B are related to connection coefficients, as we expected. 
Substituting Eq. ( |6.42|) into Eq. ( |6.38|) we have 



IQEi 



dSr 



57T 



-(E-E + B-B) 



4tt' 



-Ex B 



0(r e 



(6.41) 

(6.42) 
(6.43) 



(6.44) 



Observe that the mysterious r 2 /9 factor has disappeared, and the analogy with electromag- 
netism is improved: E and B now have their usual dimension (L" 1 ), as do the energy density 
and Poynting flux terms. I emphasize again that, in contrast to e and b in Eq. (|6.38|) , E and 
B are fields measured by observers residing in S, in the true quasilocal spirit. E is clearly 
associated with tidal forces tangential to S, and B is a measure of frame dragging effects. 
Notice that E and B vanish as r — ► 0, in accord with the equivalence principle. 

To conclude this section we consider the connection in the normal bundle and its associ- 
ated curvature. I find that 



A, 



where *B\ a 



j K B 1K , and 



\r aExj + -*B XJ +0(r 2 ), 



JF = 2£n + 0(r). 



(6.45) 
(6.46) 
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So to leading order the curvature of the normal bundle is (twice) the radial-radial component 
of the magnetic part of the Weyl tensor, and is thus associated with 'gravitational magnetic 
charge'. There are both local and global dimensions to this result. Locally, T is associated 
with frame dragging, a ready example being T for the Kerr black hole given in Eq. ( |5.15| ), 
which is proportional to the angular momentum. Globally, it is known that in exact analogy 
with the scalar curvature TZ, the integral of T over S is proportional to the Euler number of 
the normal bundle [pTfl . For a Euclidean-signature spacetime the normal bundle is an £0(2) 
(rather than S0(1, 1)) bundle, and there can be a nontrivial winding number, corresponding 
to a gravitational magnetic monopole. In the case of the Kerr spacetime there is no monopole 
present since, as is obvious from inspection of Eq. ( 5.15Q , the integral of T is zero. It might 



be interesting to explore topologically nontrivial cases in the context of the IQE. 

The result in Eq. fl6.46|) can actually be obtained immediately by inspection of Eq. (|2.16| ), 
assuming that the shear terms are higher order in r than JF is. To lowest order in r we then 
see that T = — 2i? i23- But -R0123 = C0123 is identically true, and thus we are led to Eq. (|6.46| ). 
So this equation is true whether or not matter is present, and is also independent of a. It is 
instructive to compare this result with the sectional curvature of S in vacuo: 

ooR = 2E U + 0(r). (6.47) 

This is essentially the same as the spatial infinity limit result given in Eq. ( p.l7|) , and is 
derived similarly. Comparing the previous two equations we see a striking electric/magnetic 
duality between the sectional curvature of S (electric), and the curvature of its normal 
bundle (magnetic). When matter is present, the right hand side of the equation above 
acquires an additional term, and it is precisely this term that is responsible for the 0(r 3 ) 
matter contribution seen in Eq. (|6.13|) . So the sectional curvature is the dominant term in 
the energy that encodes information about the matter content of the spacetime. It seems 
reasonable to expect inertial effects (frame dragging) produced by this matter to also play 
a role in the energy. But consideration of such effects is subtle, because the magnetic part 
of the Weyl tensor has no Newtonian gravity analogue. I have argued that the procedure 
suggested in Sec. 4 is a geometrically natural way to incorporate such inertial effects into the 
energy: one demands jF ref = J 7 , and then solves the embedding equations for the reference 
shear term, present in the reference energy. In this way the inertia information contained 
in T makes its presence felt in the energy. Moreover, by inspection of the purely spatial 
[a = 0) case of Eq. ( |6.45| ), one observes that the set of magnetic quantities: J 7 , A, and b, 
precisely encode the five independent components of the magnetic part of the Weyl tensor. 
It seems likely that the phenomenon of gravitational energy is subtle enough to be sensitive 
to this full set. Out of this set, in this section we have seen only the role of b. To see whether 
or not the other components play a role (via (k 2 — / 2 ) ref ) will have to wait until a solution 
to the embedding equations is found. 



7 Asymptotically anti-de Sitter spacetimes 

In this last section we will explore the significance of the aaR rei term in Eq. 
Suppose our physical spacetime (M, g) is asymptotically anti-de Sitter space. The aaR ref 
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term in IQE rcf gives us the freedom to specify the Riemann tensor of a reference spacetime, 
which in this case is naturally the Riemann tensor of anti-de Sitter space. Thus, according 
to Eq. (|OD we have aaR iei = -2/£ 2 , and so 



IQE 



ref 



— / dS 

8tt Js 



(7.1) 



In the large sphere limit it is clear that the cosmological constant term will dominate, rather 
than TZ, and the behavior of the IQE is qualitatively different from that for asymptotically 
flat spacetimes. 

Let us now specialize to the case that (M, g) is the AdS4-Schwarzschild spacetime, so that 
our main argument is not obscured by consideration of the shear terms, which will obviously 
be just zero. The line element in this case is given by 



ds 2 



N 2 dt 2 + —dr 2 + r 2 dtt 2 , where N(r) = f(r) = — + 1 



2M 



/ 



(7.2) 



and dfl 2 = d6 2 + sin 2 ^ d(p 2 is the line element on the unit round sphere. Let S be a t, r = 
constant two-sphere. Its scalar curvature is TZ = 2/r 2 , and a simple calculation shows that 
its sectional curvature is given by 



2 AM 
ooR = — — H — . 



(7.3) 



the dominant term coming from the anti-de Sitter 'background'. Substituting these results 
into Eqs. O and (Q we find 



IQE = -■£-[ 

87T JS 



dS 



2 

T 2 



AM 




^ + 0(1). (7.4) 



The divergent terms due to the cosmological constant cancel, so the limit of the IQE as 
r — > oo exists, and this limit is zero. This would be the expected result if the IQE had the 
interpretation of an energy, which should be red shifted to zero by the cosmological horizon. 
In contrast, we do not expect an invariant mass to be redshifted. This is why the terminology 
'invariant quasilocal energy' was chosen rather than 'invariant quasilocal mass', even though 
the IQE is the analogue of the mass m in the formula: m = y/E 2 — p 2 . 

However, one can easily modify the definition of the IQE — to give it the interpretation of 
mass — by multiplying the right hand side of Eq. (|3.4j ) by a lapse function. Thus one replaces 
Eq. O with 



(7.5) 



IQE[iV B ] = / dSN B J2[TZ-aaR} + ^- [ dSNfJ2 [TZ - aaR^} 
on Js on Js 



(ignoring the shear terms). Here the smearing function, Ng, is the lapse function in the 
timelike three-boundary, B, swept out by the two-parameter family of observers (cf. Eqs. (11) 
and (13) in Ref. |P6|| ). In the AdS4-Schwarzschild example B is an r = constant surface, and 
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Njg = N. The question arises, What are we to put for 7Vg ef ? The answer that works is 
Ng = N B , which is intuitively justified as follows: we are already isometrically embedding 
(S,a) into (M ref , g ref ), and the condition N B cf = N B represents the next would-be step 
towards an isometric embedding of (-6,7) into (M ref , g mf ), where 7^ is the three-metric in 
£>. ('Would-be' in the sense that, while the lapse carries some information about 7^, we 
still only need to embed (S, a) into (M ref , g Tef ), not (23, 7) into (M ref , g rei ).) By comparing 
Eq. (|7.5| ) with ( |7.4|) , and using the fact that the lapse function goes as r/£ for large r, it is 
easy to see that with N B cf = N B we get lim^oo IQEfiVg] = M for the AdS/i-Schwarzschild 
case. Thus IQE[iVs] has the interpretation of a mass, as claimed. Unlike the original IQE, it 
is not red shifted to zero, and is thus a different physical quantity. Regarding the comment 
at the end of the previous paragraph, since special relativity does not know about lapse 
functions, the generalization given in Eq. ( |3.4|) is open to this ambiguity: one can define 
both an invariant quasilocal energy and an 'invariant' quasilocal mass. 

It is instructive to evaluate IQE[iVg] also at the horizon, r = r + (where N B (r + ) = 0), 
and compare with what one gets using the unsmeared IQE. The following results for the 
AdSzi-Schwarzschild example are easily established: 



IQE= Y * Z * (7-6) 




IQEM = {° M % T T Z r * P.7) 

We thus learn that the IQE decreases with r, which can be interpreted as the result of 
negative binding energy — another reason to think of the IQE as an energy. On the other 
hand, IQEfiVg] increases with r, which might be interpreted as saying that, for larger r, 
'more mass is enclosed'; it starts from zero at the horizon (no mass in the interior of the 
black hole) and accumulates to M at infinity. This is reminiscent of the old notion that the 
substance of mass is nothing but the curvature of spacetime itself. A similar behavior is 
observed for the usual Schwarzschild case: 

|2Mat r = 2M 

\ M at r = 00 v ' 

T „_. r , r , f at r = 2M ,„ 
IQE[N B ] = I M atr = oQ (7.9) 

The only qualitative difference occurs at r = 00, where IQEfiVg] = IQE in the Schwarzschild 
case because, of course, the lapse function goes to one in this limit. There is no cosmological 
horizon. 

Despite these appealing features of IQE[iVg], it is unsatisfactory from the point of view 
taken here because the presence of the lapse function means it depends on a choice of three- 
surface passing through S (hence the use of inverted commas in the terminology 'invariant' 
quasilocal mass above). The situation might be improved by replacing k with N B k and / with 
N^l in Eq. ( |3.4j ), and then proceeding as before. Here N B and N s are time and radial lapse 
functions, respectively (equal to iV and 1/ f in the AdS4-Schwarzschild example) . Admittedly 
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such a procedure is ad hoc, and unless it can be improved upon we are not particularly 
interested in IQEfiVg]. It was introduced simply to illustrate the distinction between mass 
and energy, but for the remainder of this section we will return to the original definition of 
the IQE. 

The main point of this section is to draw attention to a remarkable similarity between the 
reference subtraction term given in Eq. ( |7. 1| ) , and a certain counterterm action recently sug- 
gested in the context of the conjectured AdS/CFT correspondence. We begin by observing 
that when the reference shear term vanishes, Eq. ( |7.1| ) reduces to precisely the same refer- 
ence subtraction term suggested by Lau f37j, in the context of the Brown- York quasilocal 
energy (except that Lau's expression has a lapse function present in the manner discussed 
above). However, our derivations of this expression are different. Lau employs a 'light cone 
reference' embedding of (S,a), together with a 'rest frame' assumption, Z ref = 0, to derive 
an expression for fc ref , which is then used to construct his reference subtraction term. In 
our case we get the same end result, but we get it without any recourse to a reference em- 
bedding! This is because the cosmological constant term in Eq. ( |7.1| ) comes directly from 
the explicit dependence of IQE ref on the Riemann tensor of the reference spacetime, i.e., the 
term aaR rei in Eq. ( |4.3| ). Our reference embedding of (S, a) into (M , g rei ) is required only 
to evaluate the reference shear term, (k 2 — Z 2 ) rcf . This 'higher order correction' — which I 
have argued accounts for angular momentum — is not present in Lau's reference subtraction 
term. Also, his additional 'rest frame' assumption is not required here because the IQE is 
already naturally a 'rest frame energy'. 

Let us now return to Eq. ( |1.5| ). We know that when space is non-compact the boundary 
(or quasilocal) energy-momentum tensor TJj b = — U ab / (8tt) diverges in general as B is taken 
to infinity^ To render it finite, Brown and York suggest the use of a reference subtraction 
term that involves an isometric embedding of (£>, 7) into a suitable reference spacetime. How- 
ever, like their prescription to embed (S, a) into a suitable three-dimensional reference space, 
this prescription suffers from the drawback that such a co dimension-one embedding does not 
always exist. Recently, Balasubramanian and Kraus [3~5 have proposed an alternative pro- 
cedure: Since it is always possible to add to the action a local functional of the intrinsic 
geometry of the boundary without affecting the equations of motion or the symmetries (but 
of course this alters T^ b ), their idea is to choose this functional such that its divergences 
cancel those of the original T^ fc , rendering the improved boundary energy- momentum tensor 
finite as B is taken to infinity. No recourse is made to a reference embedding. This procedure 
was first applied to spacetimes that are asymptotically AdS space, in which case the required 
'counterterms' amounted to a simple finite polynomial in the curvature invariants of B . 
This idea is exactly analogous to the standard prescription for removing ultraviolet diver- 
gences in quantum field theory by adding to the Lagrangian a finite polynomial in the fields. 
Moreover, the conjectured AdS/CFT correspondence [47| implies that the two procedures 
are not merely analogous, they are one and the same |36] . 

Now since flat spacetime is recovered from AdS space by taking i to infinity, one might 

19 As elsewhere in this paper, we use the symbol B loosely to refer to either a timelike three-surface in the 
interior of M, bounding a finite spatial region, or the boundary at infinity. It's meaning should be clear from 
the context in which it is used. 
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expect that in this same limit the counterterms found by Balasubramanian and Kraus would 
produce counterterms suitable for asymptotically flat spacetimes. It is not obvious that this 
is so |36| . However, Mann [ 48| has suggested the following generalization of their counterterm 
action: 



57T JB, 




(7.10) 



where -R(7) is the scalar curvature of (£>, 7), and B^ indicates that we are to take the limit as 
B goes to infinity. For small I Mann's formula reduces to the one given by Balasubramanian 
and Kraus, but in addition it has a smooth flat spacetime limit as i — > 00. Moreover, Mann 
showed that in many explicit examples it leads to a cancellation of all divergences, and the 



remaining finite part agrees with that obtained using the reference spacetime procedure 
49fl . While a counterterm action and a reference mass are not the same thing, the resemblance 
between the expressions in Eqs. ( 7.10 ) and (7.1) is nevertheless striking.^] 



To see that the connection between IQE re in Eq. (|7. 1|) and the AdS / CFT-inspired coun- 
terterm action is probably much deeper than mere resemblance, we now turn to recent work 
done by Kraus et al ||50|| . Besides providing an independent derivation of Mann's formula, 
and its generalization to higher dimensions, of most interest to us here is their geometrical 
argument suggesting what the counterterm for Il ab in Eq. (L5) should be in order to cancel 
divergences. Their result is an expansion in powers of £. Denoting their counterterm (fl a b) 
as ITfL and specializing their result to a three-dimensional boundary B, they find: 



n 



ct 
(ib 



~ n lab + £[Rab - -labR 



] \i ab (R cd R cd -\R 2 



+- A RR ab - 2R cd R acbd + - A D a D b R - UR ah + ~ 7a6 Di?} + 0(£ 5 ). (7.11) 

The curvatures and covariant derivatives in this expression all refer to the induced timelike 
three-metric •y ab on B. 

Now consider our usual two-surface (S, a) in the physical spacetime (M,g), which is an 
asymptotically AdS space. As we did at the end of Sec. 4, suppose that S is such that 
(k 2 — I 2 ) > and k > 0. Then we can always find a timelike unit vector u a normal to S such 
that / = 0, and so \/k 2 — I 2 = k = U ab u a u b is the Brown- York energy surface density (modulo 
the factor of — l/(87r)). In other words, our unreferenced IQE reduces to the unreferenced 
Brown- York CQE, which would be called the 'unrenormalized' energy in Ref. 
counterterm required to renormalize the energy surface density is thus IT 
will denote at E ct . Hence, our task is to compare E ct := H 



50 . The 
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ct „,a„,b 



fib 



U U 



which we 



ab 



U U 



with the integrand of IQE 



ref 



in Eq. (|7. 1|) ; we expect to see at least some measure of agreement between the two. 

This comparison will not be straightforward, however, because on the one hand we expect 
the integrand of IQE ref to depend on 1Z, JF, and their derivatives in S, as discussed previously, 
whereas on the other hand E ct depends on the three-metric 7^. Nevertheless, let us see how 
far we can get. Let B be a three-surface in (M,g) passing through S in a direction tangent 
to u a on S. Different choices of B satisfying these conditions will lead to different induced 

20 I am indebted to R.B. Mann for pointing out to me the significance of the aaR rei term in fQE rc t , an d 
emphasizing that it provides at least some measure of geometrical motivation for his expression in Eq. ( [7.10 ). 
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metrics jab, but this ambiguity will not affect our considerations. At least jab on S is uniquely 
determined, and some information about 7^ in the neighborhood of S is determined by the 
condition I = 0. Our choice of B means that l a b defined in Eqs. ( |2.8|) is the extrinsic curvature 
of (S, a) as embedded in (£>, 7), and so the corresponding co dimension-one Gauss embedding 
equation reads 



efgh 



T^abcd iJ'aJ'bd ^bJ 



bc^ad ) 



This is just a truncated version of Eq. 
(B, 7 ),not(M,(/). 



(7.12) 

except here R e f g h is the Riemann tensor of 



Now let denote the term in E ct of order £ n . Inspection of Eq. ( |7.11|) shows that 
2/1. The term can be written in terms of G a b, the Einstein tensor of (£>, 7), 



(-1) 



E, 



and we have 



E. 



(i) 



ct 



£G ab u a u b 



— CT <J K, 



abed 



2<* 



(7.13) 



The second equality is an easily derived identity (valid in any codimension-one setting) relat- 
ing the uu component of the Einstein tensor to the sectional curvature of the hypersurface 
orthogonal to u a (in this case the hypersurface S in B). The third equality follows from 
contracting Eq. (|7.12| ) with a ac <j bd , and using the fact that / = by our choice of u a . (And 
as usual, I 2 is shorthand for l a d ab -) Thus, to order £ we have 



E, 



ct 



7 + 2 (K 



P) + 0(£ 3 



P + 0(£ 2 



(7.14) 



Comparing the last expression with Eq. ( |7.1| ) suggests the correspondence: 



(P - I 



2\re{ 



-P + Oil 2 



(7.15) 



Immediately we see something odd: we are identifying a boost invariant quantity with one 
that is not, i.e., it seems that a k 2 is missing from the right hand side. I will comment 
on this shortly. Let us assume for the moment that the right hand side reads (k 2 — I 2 ), in 
which case Eq. ( |7.15| ) seems reasonable: it suggests that, if we solve the embedding equations 
( |4.6[ - [4.8|) for (k 2 — / 2 ) ref we will find that, to lowest order in £, the reference shear term is 
the same as the unreferenced shear term, the difference to be seen at a higher order in I. 
On the other hand, this seems like a problem: Would it not mean, e.g., that the shear 
terms in Eq. ( |5.5| ) basically cancel, thus ruining the Bondi-Sachs mass result in Eq.( |5.35| ), 
which depends so crucially on (k 2 — Z 2 )? The answer is No, because Eq. ( |5.5|) is valid in 
the asymptotically flat case, not the asymptotically AdS case. To make a statement that is 
valid in the asymptotically flat case (£ — * 00) we need to know E ci to all orders in £, then 
sum the infinite series, and finally take the limit I — > 00. So being at the 'other end' of the 
series, Eq. ( |7.15| ) has nothing to say about the asymptotically flat case. But we also expected 
ik 2 — l 2 ) rei to depend on 71, J 7 , and their derivatives. Why do we not see these quantities 
on the right hand side of Eq. ( |7.15| )? The answer is, We will — we just have to calculate E ct 
to the next order in £. 

But before doing so I will comment on the 'missing' k 2 in Eq. ( |7.15[ ). Kraus et al pU[ have 
devised an algorithm to compute the extrinsic geometrical quantity 11^ from the intrinsic 
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geometry of B. Insofar as l ab (and thus I 2 ) depends only on the metric jab, there is no doubt 
that the —I 2 term in Eq. (|7.15|) is 'correct'. On the other hand, k ao (and thus k 2 ) depends 
on the extrinsic geometry of B, being just a certain projection of U ao into S. The algorithm 
of Kraus et al relies on the fact |jTJ that the divergent part of the derivative of Tl ao in the 
direction normal to B can be expressed in terms of just the intrinsic geometry of B. In 
essence, their algorithm is designed precisely to compute the divergent part of U a b- The 
'correctness' of the accompanying finite part is a subtle issue. In a slightly different context, 
they discuss two different counterterm actions that both properly cancel divergences, but 
that lead to different finite terms in the action. Furthermore, they point out that their 
algorithm, when carried to all orders in £, might imply singularities in the bulk spacetime, 
but that this is of no concern because they truncate their counterterm expressions to a 
finite number of terms, enough at least to cancel the divergences. In our case we have the 
quasilocal idea in mind, and so are interested in all of the finite terms — it matters what 
happens in the bulk. But going further with this discussion will take us beyond the scope 
set for this simple comparison. I will just conclude by saying that, insofar as the shear terms 
almost certainly represent a finite contribution to the energy, we do not necessarily expect 
the algorithm of Kraus et al to produce a k 2 term on the right hand side of Eq. ( |7.15| ). Our 
goals are slightly different, and it is too much to expect exact agreement between E ct and 
the integrand of IQE ref . 

Nevertheless, it is still instructive to proceed with the comparison to the next order in 
I. In light of my previous remarks, we will make the simplifying assumption that the metric 
on B has a product structure: j a b dx a dx b = —N 2 dt 2 + Uij(x) dx l dx\ where x a = (t, x l ) are 
local coordinates on B, N is a constant lapse, and <Jij{x) is the metric on any t = constant 
two-surface 5*. The idea is that E ct and the integrand of IQE ref should agree at least in their 
dependence on the intrinsic geometry of (S,a). Assuming such a product structure for j ab 
is a convenient was to isolate this dependence, and ignore everything else. 

In this case the only nonvanishing components of the Riemann tensor of j a b are Rijki = 



the Riemann tensor of cx,- 



£lZ/2, and it is a simple exercise to work out E t 



And clearly l ab = 0. 

(3) 
ct 



Thus in Eq. ( |7.13| ) reduces to 
The net result is 



E 



ct 



+ 2^~16^ 2 " 4A ^) +0(£5) = \ 



(7.16) 



where A is the Laplacian in (S, ex). Comparing the last expression with Eq. we now 
have the higher order correspondence: 



(P - l 2 y cf < — > -— ATZ + 0(£ A ). 



(7.17) 



Thus we begin to see how a solution to our embedding equations might yield an expression 
for (k — I 2 ) in terms of 1Z, J 7 , and their derivatives, as we have expected all along. 

To conclude this section we make two general observations. First, it is especially clear 
from the higher order expression in Eq. (|7. 16| ) that the AdS/CFT-inspired counterterm 
energy is, in fact, the square root of some quantity. This is not surprising, since the algorithm 
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of Kraus et al is a means of solving a Gauss embedding equation for IL^, and this 
equation is quadratic in TV£ b . But it is significant. Beginning simply with the definition of 
the quasilocal energy-momentum tensor as the functional derivative of the action with respect 
to the boundary metric 0, in which there is no square root in sight, the counterterm energy 
required to cancel divergences unmistakably involves a square root. Moreover, it concurs 
with the square root introduced here, in the context of the IQE, as the general relativistic 
analogue of the special relativistic formula: m = y/E 2 — p 2 . I believe it is unlikely this is 
a mere coincidence. Given that it is nonanalytic, a square root is too unusual an object to 
occur without good reason. 

Secondly, under the square root (in our case) is 2/£ 2 + 1Z + (k 2 — l 2 ) re{ . In the case of 
the AdS/CFT-inspired counterterm energy ||50|| , it is 2/£ 2 + TZ + X, where X is an infinite 
series in increasing powers of £. That X is clearly not zero lends strong support for our 
additional term [k 2 — Z 2 ) rcf , which is thus seen to be a necessary generalization of Lau's 
suggestion |37|. I have argued that its necessity is closely linked to the proper inclusion 
of angular momentum in the energy. Given that angular momentum is a subtle notion in 
general relativity, especially so at the quasilocal level we envision here, it is not surprising 
that our biggest difficulty lies in evaluating (k 2 — / 2 ) ref . In light of the algorithm given by 
Kraus et al [5(J, work is currently in progress to try to apply similar techniques to solve 



the embedding equations ( |4.6| - |4.8[ ). Since these embedding equations are manifestly boost 
invariant, I expect at least to recover the 'missing' k 2 term in Eq. ( |7.15| ), and hopefully the 
entire series. 



8 Summary and discussion 

In this paper I have introduced a new definition of quasilocal energy that is a simple 
modification of the Brown- York quasilocal energy. I just replace their energy surface density 
k with \Jk 2 — I 2 , where I is the radial momentum surface density. (For ease of exposition 
here I will omit the — l/(87r) factors.) The principle motivation for doing this stems from 
an analogy with the formula: m = \JE 2 — p 2 in special relativity. Identifying E with k 
(which are both energies), and p with I (both momenta), identifies m with \Jk 2 — I 2 . Like 
m, \Jk 2 — I 2 is a boost invariant quantity, and hence the integral of \Jk 2 — I 2 over a spacelike 
two-surface S gives rise to an 'invariant quasilocal energy', or IQE. In what follows I will 
refer to the Brown- York quasilocal energy as the CQE — canonical quasilocal energy. 

There are several important consequences of replacing k with \/k 2 — I 2 : 

• While k is always well defined for any spacelike two-surface S, \/k 2 — I 2 is not. Roughly 
speaking, it is real when S lies in the exterior region of a black hole, zero when it is on 
the horizon, and imaginary in the black hole interior. Thus (again roughly speaking) 
the IQE asserts that energy is real only outside of a black hole. 

• Both the CQE and the IQE require a reference energy subtraction procedure. Since k 
is associated with a spacelike three-surface spanning S, the reference space into which 
S is to be isometrically embedded is inherently three-dimensional. Such a codimension- 
one embedding does not always exist, but when it does, it is essentially unique. This 
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means the CQE, when it is defined, is unique. In contrast, \Jk? — I 2 makes no ref- 
erence to a three-surface spanning S, and so the reference space(time) is inherently 
four-dimensional. Such co dimension-two embeddings (at least of a generic non-round 



sphere into Minkowski space) always exist p5|, but are not unique. However, in this 



situation there are two curvatures associated with S: the scalar curvature 71, and the 
curvature of the normal bundle, T . A necessary condition for an isometric embedding 
is that lZ rei = JZ. I argued that demanding also ^-" ref = T is both a means to make 
the embedding essentially unique, and at the same time, a geometrically natural way 
to properly incorporate angular momentum into energy at the quasilocal level. In- 
deed, since angular momentum is associated with rotational kinetic energy, it should 
contribute to the energy in some way. 

While CQE ref is associated with a reference energy density fc ref , IQE ref is concerned 
with a reference shear term, {k 2 — l 2 ) rei . (k a b and l a b are the trace-free parts of the 
two extrinsic curvatures of S.) In a certain sense, the IQE already inherently contains 
the correct reference energy, without recourse to a reference embedding. The reference 
embedding is required only to determine the reference shear term, which is a higher 
order correction to the energy associated with angular momentum. 

The CQE is sensitive to the sign of k, whereas since it involves \/k 2 — I 2 , the IQE 
is not. Thus one can easily construct simple examples for which the two energies 
give different results, even when / = 0. Thus the IQE is not simply the 'rest energy' 
version of the CQE. Note: the IQE naturally assigns zero energy to any two-surface 
in flat spacetime. This is because the natural reference spacetime in this case is the 
very same spacetime, namely flat spacetime. So obviously one can always reference- 
embed the two-surface identically (up to Poincare transformations) to the way it is 
embedded in the physical spacetime, and get IQE = 0. The only subtlety that may 
arise is if the two conditions: 'isometric embedding' and 'jF ref = JF' do not uniquely 
determine (k 2 — Z 2 ) ref . Then the flat spacetime result (IQE = 0) may be reduced 
to a choice, rather than a necessary fact. To properly address this subtlety requires 
an in depth understanding of the embedding equations. But in any case, the fact 
that IQE = in flat spacetime is independent of the motion of the observers. In 
contrast, moving observers in flat spacetime could measure nonzero energy in the 
Brown- York approach This is because under a radial boost the Brown- York 



energy surface density dilates by a Lorentz factor, as in special relativity, whereas the 
reference energy surface density does not. According to Ref. |7j the latter depends only 
on the intrinsic geometry of S, and therefore does not know about the time derivative 
of this geometry^ 

We examined both the large and small sphere limits of the IQE, taking S to be asymptot- 
ically round for simplicity. In an asymptotically flat spacetime, the large sphere limit of the 
IQE in a spatial direction yields the ADM mass. In the future null direction it reduces to the 
Bondi-Sachs mass, provided the reference shear term is a total divergence. Short of solving 



21 I thank R.B. Mann and I.S. Booth for this remark on the Brown- York case. 
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the embedding equations, I gave a heuristic argument which shows that is is. It is significant 
that this argument relies on the condition jF rcf = JF, since this provides evidence that the 
curvature of the normal bundle is involved in quasilocal energy, albeit its involvement in this 
simple example is minimal. 

The quantity y/k 2 — P is proportional to the mean curvature of S as a two-surface em- 
bedded in the physical spacetime, and so the IQE is a natural geometrical invariant of S. 
Since the Hawking mass [25] is constructed using (k 2 — I 2 ), the IQE can be thought of 
roughly as the square root of the Hawking mass. In the small sphere limit the square root 
disappears, and to leading order the IQE reduces to the Hawking mass (but differs from it 
at higher order). Thus, when matter is present, the lowest order contribution to the IQE 
gives the standard result: (47rr 3 /3) T^ at u a u b , i.e., the expected matter energy contained in 
a small sphere of proper radius r. Note that u a here is not necessarily the four-velocity of 
any observer on S, since the IQE is boost invariant, and so independent of the observers' 
velocities on S. Rather, u a is the four- velocity that observers would have i/they were in the 
rest frame determined by S. More precisely, in the small sphere limit we considered, namely 
a t, r — constant two-sphere in Riemann normal coordinates (with t oc r), u a = (d/dt) a 
evaluated at the center of the sphere. In the limit r — > 0, the four- velocity (d/dt) a corre- 
sponds to observers who at each point on S have zero radial momentum, i.e., 1 = 0. In 
general, since the IQE is an energy rather than a mass, the question arises, In whose rest 
frame is the energy measured?^ The answer is, The 'quasilocal rest frame' determined by 
the condition / = at each point on S. Whenever (k 2 — I 2 ) > 0, observers on S can always 
achieve this state of motion by appropriate local radial boosts. This (or more precisely, 
(k 2 — I 2 ) > 0) is the same condition required for the unreferenced IQE to be well defined in 
the first place — refer to the text in the first bullet on page ^7. 

Returning to the small sphere case, in vacuo the leading order contribution due to gravi- 
tational energy occurs at order r 5 . At this order the IQE results are inconclusive because it 
is expected that the reference shear term will play a significant role, and without a solution 
to the embedding equations (which is an extremely difficult problem) this term cannot be 
determined. Nevertheless, it was possible to show that in the small sphere limit, the Hawking 
mass, which in this case is closely related to the IQE, can be understood as a measure of 
the gravitational energy contained in S by considering certain tangential 'electromagnetic' 
fields E and B induced on S by the Weyl curvature S is embedded in. In terms of E and B, 
gravitational energy and radiation are essentially identical in nature to their counterparts 
in electromagnetism, except for one crucial difference: the density r(E ■ E + B ■ B)/(8%) is 
integrated over the surface S to determine the gravitational energy contained in the spatial 
volume that 5* encloses. Here r can be thought of as the areal radius of S, so the measurement 
is truly quasilocal. 

The IQE was analyzed in the context of asymptotically anti-de Sitter spacetimes. The 
fact that IQE ref depends explicitly on the Riemann tensor of the reference spacetime (nat- 
urally taken to be anti-de Sitter space) was seen to play a significant role. A connection 
was established between IQE ref and a certain counterterm energy that has recently been 
proposed p0| in the context of the conjectured AdS/CFT correspondence. Two similarities 



2 I thank A. Ashtekar for posing this question. 
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are striking: (i) Both energies involve a square root, and (ii) the two leading terms under the 
square root match. The remaining term under the square root in our case is the reference 



shear term, [k — I ) ; in the case of Ref. |50| it is an infinite series in increasing powers 



of l 1 where I is the radius of curvature of the AdS space. It was shown that the first two 
nontrivial terms of this series (i.e., to the highest order given in Ref. [|50|) can plausibly be 
identified with (k 2 — Z 2 ) ref . This agreement is impressive because IQE ref and the AdS/CFT- 
inspired counterterm energy are independently motivated, and derived quite differently. It 



might be possible to use techniques developed in Ref. |50| to solve our embedding equations 
for (k 2 — Z 2 ) ref . The present lack of a solution to these equations is the main outstanding 
obstacle to further understanding the nature of the IQE. 

A final remark is in order. Most definitions of quasilocal energy, including the IQE, 
assume that energy is associated with a closed spacelike two-surface, S. Given such a two- 
surface one can always find a timelike unit normal vector field u a , which at each point on 
S is supposed to correspond to an observer's instantaneous four- velocity. But this may not 
be a general enough setting. While a two-parameter family of observers will always sweep 
out a timelike three-surface B, the two-surface elements orthogonal to their world lines in 
B are not, in general, integrable. Thus a shift in emphasis from S to B, i.e., from Eulerian 



to Lorentzian observers p8| , might lead to a deeper understanding of quasilocal energy, 
in particular of gravitational radiation at the quasilocal level. This shift would also bring 
the quasilocal energy idea closer in line with the conjectured AdS/CFT correspondence. 
Whether or not this is the right direction, the results in Sec. 7 strongly suggest that this is 
the direction the IQE is pointing in. 
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